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Abstract 

Searching for heavy neutral gauge bosons Z' , predicted in extensions of the 
Standard Model based on a U(l)' gauge symmetry, is one of the challenging ob- 
jectives of the experiments carried out at the Large Hadron Collider. In this paper, 
we study Z' phenomenology at hadron colliders according to several U(l /-based 
models and in the Sequential Standard Model. In particular, possible Z' decays into 
supersymmetric particles are included, in addition to the Standard Model modes so 
far investigated. We point out the impact of the U( 1 )' group on the MSSM spectrum 
and, for a better understanding, we consider a few benchmarks points in the param- 
eter space. We account for the D-term contribution, due to the breaking of U(l)', 
to slepton and squark masses and investigate its effect on Z' decays into sfermions. 
Results on branching ratios and cross sections are presented, as a function of the 
MSSM and U(l)' parameters, which are varied within suitable ranges. We pay spe- 
cial attention to final states with leptons and missing energy and make predictions 
on the number of events with sparticle production in Z' decays, for a few values of 
integrated luminosity and centre-of-mass energy of the LHC. 
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1 Introduction 



The Standard Model (SM) of the strong and electroweak interactions has been so far success- 
fully tested at several machines, such the LEP and Tevatron accelerators and has been lately 
confirmed by the data collected by the Large Hadron Collider (LHC). New physics models 
have nonetheless been proposed to solve the drawbacks of the SM, namely the hierarchy prob- 
lem, the Dark Matter observation or the still undetected Higgs boson, responsible for the mass 
generation. The large amount of data collected at the centre-of-mass energy of 7 TeV at the 
LHC opens a window to extensively search for new physics. The further increase to 8 and ul- 
timately 14 TeV, as well as higher integrated luminosities, will extend this investigation in the 
near future. 

The simplest possible extension of the SM consists in a gauge group of larger rank involving 
the introduction of one extra U(l)' factor, inspired by Grand Unification Theories (GUTs), 
which leads to the prediction of a new neutral gauge boson Z' . The phenomenology of the Z' 
has been studied from a theoretical viewpoint (see, e.g., the reviews [HJI3 or the more recent 
work in Refs. (3M1), whereas searches for new heavy gauge bosons have been carried out at the 
Tevatron by the CDF and DO HI Collaborations and at the LHC by ATLAS and CMS M- 
Besides the Z' bosons yielded by the extra U(l)' group, the analyses have also investigated the 
so-called Sequential Standard Model (Z' SSM ), i.e. a Z' with the same couplings to fermions and 
gauge bosons as the Z of the SM. The Sequential Standard Model does not have theoretical 
bases like the U(l)' models, but it is used as a benchmark, since, as will be seen later on, the 
production cross section is just function of the Z' mass and there is no dependence on other 
parameters. 

The Tevatron analyses searched for high-mass dielectron resonances in pp collisions at 1 .96 
TeV and set a lower Z' mass limit of about 1023 (DO) and 963 (CDF) GeV for the Z' SSM . The 
LHC experiments investigated the production of both dielectrons and dimuons at large invariant 
masses and several models of Z' production, i.e. different U(l)' gauge groups. The CMS 
Collaboration, by using event samples corresponding to an integrated luminosity of 1.1 fb _1 , 
excluded a Z' with SM-like couplings and mass below 2.32 TeV, a GUT-inspired Z' below 1.49- 
1.69 TeV and a Kaluza-Klein graviton in extra-dimension models [9] below 0.71-1.63 TeV. 
The ATLAS Collaboration analyzed 5 fb 1 of data and obtained a bound of 2.21 TeV for the 
SM-like case, in the range 1.76-1.96 TeV for the U(l)' scenarios and about 0.91-2.16 TeV for 
the Randall-Sundrum gravitons 0. 

All such analyses, and therefore the obtained exclusion limits, crucially rely on the assump- 
tion that the Z' decays into Standard Model particles, with branching ratios depending on its 
mass and, in the GUT-driven case, on the parameters characterizing the specific U(l)' model: 
such a choice is dictated by the sake of minimizing the parameters ruling the Z' phenomenology. 
As a matter of fact, in the perspective of searching for new physics at the LHC, there is no actual 
reason to exclude Z' decays into channels beyond the SM, such as its supersymmetry. In fact, 
new physics contributions to the Z' width will significantly decrease the branching ratios into 
SM particles, and therefore the mass limits quoted by the experiments may have to be revisited. 
Furthermore, Z' decays into supersymmetric particles, if existing, represent an excellent tool 
to investigate the electroweak interactions at the LHC in a phase- space comer that cannot be 
explored by employing the usual techniques. Therefore, the possible discovery of supersym- 

1 *The exclusion ranges depend on the specific U( 1 )' model and, for the graviton searches, on the coupling value. 
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metry in Z ! processes would open the road to additional investigations, since one would need to 
formulate a scenario accommodating both sparticles and heavy gauge bosons. 

The scope of this paper is indeed the investigation of the phenomenology of Z 1 bosons at 
the LHC, assuming that they can decay into both SM and supersymmetric particles. As for 
supersymmetry, we shall refer to the Minimal Supersymmetric Standard Model (MSSM) ifTOl 
[TT1| and study the dependence on the MSSM parameters. A pioneering study of supersymmetric 
contributions to Z' decays was carried out in [fT2j. wherein the partial widths in all SM and 
MSSM channels were derived analytically, and the branching ratios computed for a few U(l)' 
scenarios. However, the numerical analysis was performed for a mass mz'=l 00 GeV, presently 
ruled out by the late experimental measurements, and only for one point of the supersymmetric 
phase space. Therefore, no firm conclusion could be drawn about the feasibility to search for 
the Z' within supersymmetry at the LHC. This issue was tackled again more recently. Ref. IfPJI 
studied how the Z' mass exclusion limits change once sparticle and exotic decay modes are 
included, for many U(l)' models and varying the supersymmetric particle masses from to 
2.5 TeV. The Higgs and neutralino sectors in extensions of the MSSM, including GUT-inspired 
U(l)' models, were thoroughly debated in lfl4l and [031 . respectively. Ref. lfT6l considered 
the U(1)'b-l gauge group, B and L being the baryon and lepton numbers, and focused on the 
decay of the Z' into charged- slep ton pairs for a few points in the MSSM phase space and various 
values of Z' and slepton masses. Ref. ifTTl investigated all possible decays of the Z' in the SM 
and MSSM, and several U( 1)' models, for two sets of supersymmetric parameters and a Z ! mass 
in the 1-2 TeV range. 

In the following, we shall extend the above work in several aspects. Special attention will be 
paid to the MSSM spectrum after the addition of the U(l)' gauge symmetry. Squark and slepton 
masses will be parametrized as the sum of a soft mass and the so-called D- and F-terms [fT8l . 
In particular, accounting for the D-term has an impact on the sfermion masses, which get an 
extra contribution driven by the U(l)' group. Higgs, chargino and neutralino masses will be de- 
termined by diagonalizing the corresponding mass matrices. A detailed study will be thereafter 
undertaken by allowing the U(l)' and MSSM parameters to run within suitable ranges, taking 
into account the recent experimental limits. Throughout this work, particular care will be taken 
about the decay of the Z' into slepton pairs, i.e. charged sleptons or sneutrinos, eventually lead- 
ing to final states with four charged leptons or two charged leptons and missing energy, due 
to neutralinos. In fact, in the complex hadronic environment of the LHC, leptonic final states 
are the best channels to perform precise measurements and searches. Slepton production in Z' 
decays has the advantage that the Z' mass is a further kinematical constrain on the invariant 
mass of the slepton pair. Moreover, the extension of the MSSM by means of the U(l)' gauge 
group provides also an interesting scenario to study Dark Matter candidates, such as neutrali- 
nos |fT9ll20ll or right-handed sneutrinos [21 j, whose annihilation or scattering processes may 
proceed through the coupling with a Z' boson. 

We shall present results for the Z' production cross sections and the branching ratios into 
both Standard Model and supersymmetric final states, thoroughly scanning the 11(1)' and MSSM 
parameter spaces, which will enable one to estimate the LHC event rates with sparticle produc- 
tion in Z' decays. We point out that, in order to draw a statement on the feasibility of the LHC 
to search for supersymmetry in Z' decays, one should also account for the Standard Model 
backgrounds. However, in assessing whether the signal can be separated from the background, 
one would need to consider exclusive final states, wherein acceptance cuts on final-state jets, 
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leptons and possible missing energy, as well as detector effects, are expected to play a role. The 
framework of a Monte Carlo generator [12211231 , wherein both signal and background events 
are provided with parton showers, hadronization, underlying event and detector simulations, is 
therefore the ideal one to carry out such a comparison. We shall thus defer a detailed investi- 
gation of the backgrounds to a future study, after the implementation of our modelling for Z' 
production and decay in a Monte Carlo code. 

The paper is organized as follows. In Section 2 we shall briefly discuss the U(l)' gauge 
group yielding the Z' boson and the particle content of the MSSM. Section 3 will be devoted to 
summarize the new features of the MSSM, once it is used in conjunction with the U(l)' group. 
In Section 4, as a case study, we will choose a specific point of the MSSM/U(1)' parameter 
space, named 'Representative Point', and discuss the MSSM spectrum in this scenario. In 
Section 5, we shall present the Z ! branching ratios into SM and BSM particles for several U(l)' 
models and in the Sequential Standard Model. We will first investigate the decay rates in a 
particular 'Reference Point' of the parameter space and then vary the U(l)' mixing angle and 
the MSSM parameters. Particular attention will be devoted to the decays into sleptons and to 
the dependence of the branching fractions on the slepton mass. In Section 6 the leading-order 
cross section for Z' production in the U(l)' scenarios and in the Sequential Standard Model 
will be calculated. Besides, the number of events with sparticle production in Z' decays will be 
computed for a few energy and luminosity phases of the LHC. In Section 7 we shall summarize 
the main results of our study and make some final remarks on the future developments of the 
analysis here presented. In Appendix A the main formulas used to calculate the Z' branching 
ratios will be presented. 

2 Modelling 2! production and decay 

As discussed in the Introduction, we shall consider extensions of the Standard Model leading 
to Z' bosons, which will be allowed to decay into both SM and supersymmetric particles. For 
the sake of simplicity and minimizing the dependence of our analysis on unknown parameters, 
we shall refer to the MSSM. In this section we wish to briefly review the main aspects of the 
models used for Z' production and decay. 

2.1 U ( 1 y models and charges 

There are several possible extensions of the SM that can be achieved by adding an extra U(l)' 
gauge group, typical of string-inspired GUTs (see, e.g., Refs. Q3I21 for a review): each model 
is characterized by the coupling constants, the breaking scale of U(l)' and the scalar particle 
responsible for its breaking, the quantum numbers of fermions and bosons according to U(l)'. 
Throughout our work, we shall focus on the U(l)' models explored by the experimental collab- 
orations. 

Among the U(l)' gauge models, special care has been taken about those coming from a 
Grand Unification gauge group Eg, having rank 6, which breaks according to: 

E 6 ^SO(10)xU(l) v , (1) 

followed by 

SO(10) ^SU(5) xU(l)' (2) 
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The neutral vector bosons associated with the U(l)'^ and U(l)^ groups are called Z' v and Z' x , 
respectively. Any other model is characterized by an angle 6 and leads to a Z' boson which can 
be expressed as q|: 

Z'(8) =Z / v cosO-Z , x smO. (3) 

The orthogonal combination to Eq. @ is supposed to be relevant only at the Planck scale and 
can therefore be neglected even at LHC energies. Another model, named U(1)L, is inherited 
by the direct breaking of E 6 to the Standard Model (SM) group, i.e. SU(2)lx U(1)y, as in 
superstring-inspired models: 

E 6 ^SMxU(l);. (4) 

The yielded gauge boson is called and corresponds to a mixing angle = arccos \/5/8 in 
Eq. ©. The model orthogonal to U(l)^, i.e. = arccos \/5/8 — k/2, leads to a neutral boson 
which will be referred to as Z[. Furthermore, in the so-called secluded model, a U(l)'s model 
extends the MSSM with a singlet field S |]24|. The connection with the Eg groups is achieved 
assuming a mixing angle = arctan(vT5/9) —n/2 and a gauge boson Z' s . 

In the Grand Unification group E6 the matter superfields are included in the fundamental 
representation of dimension 27: 

27 = (<2,« c ,/,L,J c ,v c ,//, J D c ,// c ,D,5 c ) L . (5) 

In Eq. ©, Q is a doublet containing the left-handed quarks, i.e. 

Q=("A, (6) 



whereas L includes the left-handed leptons: 



(7) 



In Eqs. © and (Q, u, d and e denote generic quark and lepton flavours. Likewise, u c L , d c L , e c L 
and v£ are singlets, which are conjugate to the left-handed fields and thus correspond to right- 
handed quarks and leptons El In the case of supersymmetric extensions of the Standard Model, 
such as the MSSM, Q, L, u c L , d c L , e c L and v£ will be superfields containing also left-handed 
sfermions. Furthermore, in Eq. (0, H and H c are colour- singlet, electroweak doublets which 
can be interpreted as Higgs pairs: 

H=m,H<=(Kl (8) 



In the MSSM, H and H c are superfields containing also the supersymmetric partners of the 
Higgs bosons, i.e. the fermionic higgsinos. Another possible description of the H and H c fields 
in the representation 27 is that they consist of left-handed exotic leptons (sleptons) Af and L, with 



2 )ln Eq. ([3]) we followed the notation in fTZl and we shall stick to it throughout this paper. One can easily 
recover the notation used in [ 1 1 by replacing — > 6 — 7C/2. 

^Following [18], the conjugate fields are related to the right-handed ones via relations like u c L = w R , 
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the same SM quantum numbers as the Higgs fields in Eq. ([8]) lfT2l Moreover, in Eq. ©, D and 
D c are exotic vector-like quarks (squarks) and S c is a SM singletj^]. In our phenomenological 
analysis, as well as in those performed in Refs. H12l|16|[T7l , leptons and quarks contained in the 
H and D fields are neglected and assumed to be too heavy to contribute to Z' phenomenology. 
We are nevertheless aware that this is a quite strong assumption and that in forthcoming BSM 
investigations one may well assume that such exotics leptons and quarks (sleptons and squarks) 
are lighter than the Z' and therefore they can contribute to Z' decays. 

When E6 breaks according to Eqs. (H) and ©, the fields in Eq. (|5]) are reorganized according 
to SO(10) and SU(5). The SU(5) representations are the following: 

10 = (Q,u c \e c ) , 5 = (L, d c ) , 1 = (v c ) , 5 = (H,D C ) , 5 = (H C : D) , 1 = (5 C ). (9) 

From the point of view of SO(10), the assignment of the fields in the representations 16, 10 
and 1 is not uniquely determined. In particular, there is no actual reason to decide which 5 
representation should be included in 16 rather than in 10. The usual assignment consists in 
having in the representation 16 the SM fermions and in the 10 the exotics: 

16= (Q,u c ,e c ,L,d c ,v c ) , 10 = (H,D C ,H C ,D) , 1 = (S C ). (10) 

An alternative description is instead achieved by including H and D c in the 16, with L and d c 
in the 10; this 'unconventional' E6 scenario has been intensively studied in Refs. Il25l - l27l and 
leads to a different Z' phenomenology. In our paper, we shall assume the 'conventional' SO(10) 
representations, as in Eq. (TTOT ). Nevertheless, it can be shown If27l that, given a mixing angle 6, 
the unconventional Eg scenario can be recovered by applying the transformation: 



9 ->• 0+arctanV15. (11) 

In fact, in our phenomenological analysis, we shall also consider the U(1)'n model leading to 
the so-called boson, with a mixing angle = arctan vT5 — k/2. According to Eq. (fTTT) . the 
Z' N model corresponds to the Z' x one, but in the unconventional E6 scenario. Tabled] summarizes 
the U(l /-based models which will be investigated throughout this paper, along with the values 
of the mixing angle 6 . 

The U(l)' charges of the fields in Eq. ©, assuming that they are organized in the SO(10) 
representations as in (flOT ), are listed in Table |2] Under a generic U(l)' rotation, the charge of a 
field <f> is the following combination of the U(l)^ and U(l)' ¥ charges: 

Q' '(®) = Q' V (&) cos 0-Q' x (<t>) sin0. (12) 

Besides the U(l)' gauge groups, another model which is experimentally investigated is the 
so-called Sequential Standard Model (SSM), yielding a gauge boson Z' SSM , heavier than the Z 
boson, but with the same couplings to fermions and gauge bosons as in the SM. As discussed in 
the Introduction, although the SSM is not based on strong theoretical arguments, studying the 
Zg SM phenomenology is very useful, since it depends only on one parameter, the Z' mass, and 
therefore it can set a benchmark for the U(l)'-based analyses. 



4 ) ln the assumption that H and H c contain exotic leptons, it is: H = yjgj anc ' = {^c 

5 ) A variety of notation is in use in the literature to denote the exotic fields in the 27 representation. For example, 



in 1 2 25 , 26] the exotic quarks D and D c are called h and h c 
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Table 1: Z ! models along with the corresponding mixing angle, as given in Eq. ©. 



Model 


e 


Z' 





Z' 


-7t/2 


Z' 


arccos \/5/8 


7' 


arctan(v / 15/9) 


7' 


arccos a/5/8 — tt/2 


7' 

N 


arctan a/15 — 7t/2 



Table 2: U(l)' charges of the fields in the representation 27 of the Grand Unification group Eg. 









g 


-1 






-1 




d c 


3 




L 


3 




e c 


-1 






-5 




H 


-2 


-2 


H c 


2 


-2 


S c 





4 


D 


2 


-2 


D c 


-2 


-2 



In the following, the coupling constants of U(1) Y , SU(2)l and U(l)' will be named gi, g2 
and g', respectively, with gi = g2tanQw, &w being the Weinberg angle. We shall also assume, 
as occurs in Eg -inspired models, a proportionality relation between the two U(l) couplings, as 
originally proposed in [|2~8~1 : 




Before closing this subsection, we wish to stress that, in general, the electroweak-interaction 
eigenstates Z and Z' mix to yield the mass eigenstates, usually labelled as Z\ and Zi. Ref. 11291 
addressed this issue by using precise electroweak data from several experiments and concluded 
that the mixing angle 9 zz i is very small for any Z' model, namely sin0 zz / ~ 10 3 -10 4 . Like- 
wise, even the ZZ' mixing associated with the extra kinetic terms due to the two U(l) groups is 
small and can be neglected If30ll . 

2.2 Particle content of the Minimal Supersymmetric Standard Model 

The Minimal Supersymmetric Standard Model (MSSM) is the most investigated scenario for 
supersymmetry, as it presents a limited set of new parameters and particle content with respect 
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to the Standard Model. Above all, the MSSM contains the supersymmetric partners of the 
SM particles: scalar sfermions, such as sleptons and (£ = e,/i,T) and squarks q (q = 
u,d,c,sj,b), and fermionic gauginos, i.e. g, W ± , Z and f. It exhibits two Higgs doublets, 
which, after giving mass to W and Z bosons, lead to five scalar degrees of freedom, usually 
parametrized in terms of two CP-even neutral scalars, h and H, with h lighter than H, one 
CP-odd neutral pseudoscalar A, and a pair of charged Higgs bosons H ± . Each Higgs has a 
supersymmetric fermionic partner, named higgsino. The light scalar Higgs, i.e. h, roughly 
corresponds to the SM Higgs. 

The weak gauginos mix with the higgsinos to form the mass eigenstates: two pairs of 
charginos (X\ and X% ) and four neutralinos (X\ , %2, X^and Z4), where %\ is the lightest and X4 
the heaviest. Particle masses and couplings in the MSSM are determined after diagonalizing the 
relevant mass matrices. Hereafter, we assume the conservation of R-parity, with the values R p 
= +1 for SM particles and R p = -1 for their supersymmetric partners. This implies the existence 
of a stable Lightest Supersymmetric Particle (LSP), present in any supersymmetric decay chain. 
The lightest neutralino, i.e. X\, is often assumed to be the LSP. 

As for the Higgs sector, besides the two Higgs doublets of the MSSM, the extra Z' requires 
another singlet Higgs to break the U(l)' symmetry and give mass to the Z' itself. Moreover, two 
extra neutralinos are necessary, since one has a new neutral gaugino, i.e. the supersymmetric 
partner of the Z', and a further higgsino, associated with the above extra Higgs. As for the 
sfermions, squark and slepton masses will get an an extra contribution to the so-called D-term, 
depending on the U(l)' sfermion charges and Higgs vacuum expectation values. As will be 
discussed below, such D-terms, when summed to the soft masses and to the F-terms, will have 
a crucial impact on sfermion spectra and, whenever large and negative, they may even lead to 
discarding some MSSM/U(1)' scenarios. 



3 Extending the MSSM with the extra U ( 1 ) ' group 

In our modelling of Z' production and decay into SM as well as supersymmetric particles, 
the phenomenological analysis in Ref. [fT2l will be further expanded and generalized. In this 
section we summarize a few relevant points which are important for our discussion, referring to 
the work in [12] for more details. 



3.1 Higgs bosons in the MSSM and U ( 1 ) ' models 

The two Higgs doublets predicted by the MSSM (<I>i and <t>2) can be identified with the scalar 
components of the superfields H and H c in Eq. ©, whereas the extra Higgs ($3), necessary 
to break the U(l)' symmetry and give mass to the Z ; , is associated with the scalar part of the 
singlet S c . The three Higgs bosons are thus two weak-isospin doublets and one singlet: 

The vacuum expectation values of the neutral Higgs bosons are given by (§f) = v\j y2, with 
vi < V2 < V3. From the Higgs vacuum expectation values, one obtains the MSSM parameter 
tan/3, i.e. 

tanj3=v 2 /vi. (14) 
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Hereafter, we shall denote the Higgs charges according to the U(l)' symmetry as: 



Q[=Q'(H) , Q' 2 = Q'(H C ) , Q! 3 = Q!{S C ). (15) 

Their values can be obtained using the numbers in Table |2] and Eq. (PT21) . 

The MSSM superpotential contains a Higgs coupling term giving rise to the well-known ji 
parameter; because of the extra field 4>3, our model presents the additional contribution W = 
Acj»icj>2<I>3, leading to a trilinear scalar potential for the neutral Higgs bosons 

V X =XA X $$$. (16) 

The parameter X in Eq. (fT6l) is related to the usual \l term by means of the following relation, 
involving the vacuum expectation value of (j)® [TT2l : 

After symmetry breaking and giving mass to W, Z and Z' bosons, one is left with two charged 
(if ± ), and four neutral Higgs bosons, i.e. one pseudoscalar A and three scalars h, H and H'^cA 
Following OTTl . the charged-Higgs mass is obtained by diagonalizing the mass mixing matrix 

y/ 2 ^l_({g 2 2 /2-X 2 )v 2 + XA x v m /v 2 {g 2 /2-X 2 ) Vl v 2 + XA x v, \ 
^ H± 2 V (g 2 2 /2-X 2 ) Vl v 2 + XA x v 3 ( g 2/2-A 2 )vi + AA A v 2 v 3 /vJ U8j 

and is given by 

*k = ^+(l-2^W (19) 
sin2/3 V 82/ 

We refer to [TT2l for the mass matrix of the CP-even neutral Higgs bosons: the mass eigen- 
values are to be evaluated numerically and cannot be expressed in closed analytical form. One 
can nonetheless anticipate that the mass of the heaviest H' is typically about the Z' mass, and 
therefore the Z' cannot decay into channels containing H' . 

The mass of the pseudoscalar Higgs A is obtained after diagonalizing its 2 x 2 mass matrix 
and can be computed analytically, as done in OTTl : 

4 = ^(l + ^sin 2 2/3Y (20) 



sin 2/3 V 4v? 



3 



where v = Jv\ + v|Zl 



point out that in [ 12] the three neutral Higgs bosons are denoted by H®, with i = 1,2,3 and the pseudoscalar 
one by P°. 

7 )ln Eqs. ( Tl9l l and ( f20l > we have fixed the typing mistakes contained in Ref. fl2l . wherein the expressions for the 
masses of charged and pseudoscalar Higgs bosons contain extra factors of 2. 
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3.2 Neutralinos and charginos 



Besides the four neutralinos of the MSSM, Xi,...,%4, two extra neutralinos are required, 
namely X5 and X®, associated with the Z' and with the new neutral Higgs breaking U(l)'. 
The 6x6 neutralino mass matrix is typically written in the basis of the supersymmetric neutral 
bosons (— iB, —iW2,—iB',^\,^2^3)- It depends on the Higgs vacuum expectation values, on 
the soft masses of the gauginos B, W3 and B', named M\, M2 and M' hereafter, and on the Higgs 
U(l)' charges Q[, Q' 2 and Q' 3 . It reads: 



.0 = 



/ 


Mi 








-jgm 


\g\V2 
2^2V2 







M 2 





\g2V\ 










M 




Q' 2 g'v2 






jgivi 


Q\g'vi 











-\g2v2 


Q' 2 g'v2 


T2 Xv3 





\ 








Q' 3 g' v 3 


T2 Xv2 


T2 Xn 








\ 



T2'^ 





(21) 



/ 



The neutralino mass eigenstates (X\ , . . . and their masses are obtained numerically after 
diagonalizing the above matrix. Approximate analytic expressions for the neutrino masses, 
valid whenever M\, M2, M',v\ and V2 are much smaller than V3, can be found in [|32l . 

Since the new Z' and Higgs bosons are neutral, the chargino sector of the MSSM remains 
unchanged even after adding the extra U(l)' group. The chargino mass matrix is given by iflOl 



/ M2 y/lmw sin /3 

yv^m^cos^ — \x 



and its eigenvalues are 



%l>%2 2 



with 



\M 2 \ 2 + \n\ 2 + 2m^T J&x 
A% = ( \M 2 \ 2 + \ll | 2 + 2ml) 2 ~ 4 lAiM 2 - m 2 w sin2/3 1 2 . 



(22) 

(23) 
(24) 



3.3 Sfermions 

In principle, for the sake of a reliable determination of the sfermion masses, one would need 
to carry out a full investigation within models for supersymmetry breaking, such as gauge-, 
gravity- or anomaly-mediated mechanisms. Studying supers ymmetry-breaking scenarios goes 
nevertheless beyond the scopes of the present work. We just point out that supersymmetry 
can be spontaneously broken if the so-called D-term and/or the F-term in the MSSM scalar 
potential have non-zero vacuum expectation values, which can be achieved by means of the 
Fayet-Iliopoulos [|33l or O'Raifeartaigh 11341 mechanisms, respectively. 

The sfermion squared masses can thus be expressed as the sum of a soft term m 2 ,, often set to 
the same value for both squarks and sleptons at a given scale, and the corrections due to D- and 

8 ^The scalar potential is given, in terms of D- and F-terms, by V((j),(j)*) = F* l F, ■ + D"D a /2, with D" = 
—g„(<p*T a <p) and F; = 8W/8<j>j, where W is the superpotential, (pi are the scalar (Higgs) fields, g a and T" the 
coupling constant and the generators of the gauge group of the theory. 
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F-terms [fT8l . The F-terms are proportional to the SM fermion masses and therefore they are 
mostly relevant for the stop quarks. The D-term can be, in principle, important for both light 
and heavy sfermions and, for the purpose our study, it consists of two contributions. A first 
term is a correction due to the hyperfine splitting driven by the electro weak symmetry breaking, 
already present in the MSSM. For a fermion a of weak isospin T^ a , weak hypercharge Y a and 
electric charge Q a , this contribution to the D-term reads: 

Am« = (T^agi-YaglXvj-vl) = (r 3 , fl -e fl sin 2 lv )mlcos2/3. (25) 

A second contribution is due to possible extensions of the MSSM, such as our U(l)' group, and 
is related to the Higgs bosons which break the new symmetry: 

Ami 2 = 8 ^Q' a {Q / A + + (26) 

where Q[, Q' 2 and Q' 3 are the Higgs U(l)' charges defined in Eq. ([TBI and @ a the charge of 
sfermion a. When dealing with the Sequential Standard Model Z ; , only the first contribution to 
the D-term, Eq. (1251) . must be evaluated. 

Left- and right-handed sfermions mix and therefore, in order to obtain the mass eigenstates, 
one needs to diagonalize the following squared mass matrix: 



' hM f LR ) 2 (M^) 2 



(27) 



The value of the soft masses and the scale at which they are evaluated are in principle arbitrary, 
as long as the physical sfermion masses, obtained after diagonalizing the matrix (|T7T ). fall within 
the current experimental limits for slepton and squark searches. Following Ref. [fT2l . we assume 
a common soft mass of the order of few TeV for all the sfermions at the Z' scale and add to it 
the D- and F-term contributions. Another possibility would be, as done e.g. in [[T8l , fixing the 
soft mass at a high ultraviolet scale, such as the Planck mass, and then evolving it down to the 
typical energy of the process, by means of renormalization group equations. 

As an example, we present the expression for the matrix elements in the case of an up-type 
squark: 

(M a LL ) 2 = ( m 0j 2 + m 2+Q_^J m | cos2j3+Aro£ (28) 



1 2 

2~ 3 



[M a RR ) 2 = ( m 0j2 + m 2 + \| cos2j8+ Amg (29) 



(M U LR Y = m u [A u — [i cot/3) . (30) 

where x w = sin 2 0w, m a LR is me ^l,r soft mass at the Z' energy scale and Af = m u A u is the 
coupling constant entering in the Higgs-sfermion interaction term. 

The dependence on m z > and on the mixing angle is embedded in the Am- 2 term; analo- 
gous expressions hold for down squarks and sleptons [12J. In the following, the up-squark mass 
eigenstates will be named as u\ and U2 and their masses as m^, and mu 2 . Likewise, d\j_, and 
Vi 2 will be the mass eigenstates for down-type squarks, charged sleptons and sneutrinos and 
their masses will be denoted by , 2 and 2 , respectively. 
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The terms in Eqs. (1281) and (1291 ), as well as the the mixing term (|30l ) are inherited by the 
F-terms in the scalar potential. As the mass of SM light quarks and leptons is very small, such 
terms are typically irrelevant and the mass matrix of sleptons and light squarks is roughly diag- 
onal. On the contrary, the mixing term Mir can be relevant for top squarks, and therefore the 
stop mass eigenstates t\ 2 can in principle be different from the weak eigenstates ?l,r- However, 
we can anticipate that, as will be seen later on, for a Z' boson with a mass of the order of a few 
TeV, much higher than the top-quark mass, even the stop mixing term will be negligible. 



4 Representative Point 

The investigation on Z' production and decays into SM and BSM particles depends on several 
parameters, such as the Z' or supersymmetric particle masses; the experimental searches for 
physics beyond the Standard Model set exclusion limits on such quantities 11351 . 

In the following, we shall first consider a specific configuration of the parameter space, 
which we call 'Representative Point', to study the Z' phenomenology in a scenario yielding 
non-zero branching ratios in the more relevant decay channels. Then, each parameter will be 
varied individually, in order to investigate its relevance on the physical quantities. 

The set of parameters chosen is the following: 

m z i = 3 TeV , 9 = arccos — ^, 

\x = 200 GeV , tan/3 = 20 , A q = A e = A x = A f = 500 GeV , 

m l = m l = m l = m l = < = < = 2-5 TeV, 

Mi = 100 GeV , M 2 = 200 GeV , M' = 1 TeV, (31) 

where the value of 6 corresponds to the Z\ model and by q and t we have denoted any possible 
quark and lepton flavour, respectively. In Eq. (|3TT ) the gaugino masses M\ and M2 satisfy, within 
very good accuracy, the following relation, inspired by Grand Unification Theories: 

jj^ftrfftr. (32) 



4.1 Sfermion masses 

The sfermion masses are given by the sum of a common soft mass, which we have set to the 
same values for all squarks and sleptons at the Z' scale, as in Eq. (ISTI) . and the F- and D-terms, 
given in Eqs. (127|)-(|30"I). The D-term, and then the sfermion squared masses, is expected to 
depend strongly on the U(l)' and MSSM parameters, and can possibly be negative and large, 
up to the point of leading to an unphysical (imaginary) sfermion mass. The F-term, being 
proportional to the lepton/quark masses, is significant only for top squarks. In Fig. Q]we study 
the dependence of squark (left) and slepton (right) masses on the U(l)' mixing angle 6. The 
symbols U\ 7 2 h,2 and stand for generic up-, down-type squarks, charged sleptons and 
sneutrinos, respectively. With the parametrization in Eq. (f3Tb , in particular the fact that the 
Z' mass has been fixed to 3 TeV, a value much higher than SM quark and lepton masses, the 
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sfermion masses do not depend on the squark or slepton flavour. In this case, even the stop 
mixing term is negligible, so that the tip masses are roughly equal to those of the other up-type 
squarks. 

In Fig.[T]the mass spectra are presented in the range — 1.2 < < 0.8: in fact, for < — 1.2 
and > 0.8 the squared masses of d.2 and V2 become negative and thus unphysical, respectively, 
due to a D-tem which is negative and large. This implies that the model Z' x , corresponding to 
= — n/2, cannot be investigated within supersymmetry for the scenario in Eq. (|3TT ). as it 
does not yield a meaningful sfermion spectrum. In the following, we shall still investigate the 
phenomenology of the Z' x in a generic Two Higgs Doublet Model, but the sfermion decay modes 
will not contribute to its decay width. From Fig. [T](left) one can learn that the masses of ui and 
di are degenerate and vary from about 2.2 to 3 TeV for increasing values of 0, whereas the «2 
mass decreases from 2.7 to about 2 TeV. A stronger dependence on is exhibited by it 
is almost zero for ~ —1.2 and about 3 TeV for ~ 0.8. The slepton masses, as shown in 
Fig. Q] (right), decrease as increases: the mass of £\ is degenerate with V\ and shows a larger 
variation (from 3.7 to 2.2 TeV) than £2 (from 2.7 to 2.2 TeV). Sneutrinos V2 exhibit a remarkable 
dependence: m.y 2 can be as high as 4 TeV for ~ —1.2 and almost zero for ~ 0.8. 




Figure 1: Dependence on the U(l)' mixing angle of squark (left) and slepton (right) masses. 

The D-term correction, and therefore the sfermion masses, is also function of the Z' mass: 
this dependence is studied for the Z[ model and the parameters set as in the Representative 
Point, in the range 1 TeV < m z i < 3.5 TeV. In Fig. [2] the squark and slepton masses are plotted 
with respect to m z /, obtaining quite cumbersome results. The masses of Si 2, di, £2 and V2 
are independent of m' z ; on the contrary, and mj are degenerate and increase from 2.5 
TeV (m z i = 1 TeV) to about 3.5 TeV (m Z i = 3.5 TeV). The mass of d.2 is mj ~ 2.4 TeV for 
m z i = 1 TeV and m j ~ for m z i = 3.5 TeV; due to the large negative D-term for d~2 squarks, 
no physical solution for mj 2 is allowed above m z / = 3.5 TeV. 

The dependence of the sfermion masses on the initial values m~ and m~, set at the Z' mass 
scale, and varied from 400 GeV to 4 TeV, is presented in Fig. |3] As expected, given Eqs. (1281) - 
(l30l) . all sfermion masses are monotonically increasing function of nP~, in the case of Si, u%, di 

and £2, being the D-term negligible, they are degenerate and approximately equal to m~ _ in the 

whole explored range. The mass of the squark J2 is instead physical only for m Q ~ > 2. 1 TeV and 
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Figure 3: Sfermion masses as a function of the initial values nP~ and nv~, set at the Z' mass scale. 
Left: squarks. Right: sleptons. 



increases up to the value ~ 3.3 TeV for m~ — 4 TeV. The masses m^, and my 2 are also 
degenerate and vary from about 2.1 TeV (m9 = 400 GeV) to 4.5 TeV (mj = 4 TeV). 

We also studied the variation of the sfermion masses with respect to tan/3, in the range 
1 .5 < tan /3 < 5, and on the trilinear coupling A f, for 1 TeV < Af < 4 TeV, but found very little 
dependence on such parameters. Moreover, there is no dependence on Mi, M2 and M', which 
do not enter in the expressions of the sfermion masses, even after the D-term correction. 

4.2 Neutralino masses 

We wish to study the dependence of the neutralino masses on the parameters playing a role in 
our analysis: unlike the sfermion masses, they depend also on the gaugino masses Mi, M2 and 
M'. Table [3] reports the six neutralino masses for the parametrization in Eq. (f3TT) . For m Z ' = 
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3 TeV, Z ! decays into channels containing the heaviest neutralino x® are not permitted because 
of phase-space restrictions, and therefore they can be discarded in the Representative Point 
scenario. Being m^o — 2.54 TeV, decays into states containing x® are kinematically allowed, 
but one can already foresee very small branching ratios. 

Table 3: Neutralino masses for a Z' mass of 3 TeV and the parameters of the MSSM and U(l)' 
set as in Eq. (|3"TI) . 



ntifO 


AT 




A4 


AS 


Afi 


94.6 GeV 


156.6 GeV 


212.2 GeV 


261.0 GeV 


2541.0 GeV 


3541.0 GeV 



Figure 0] presents the dependence of the mass of the four lightest neutralinos, i.e. X\ ■ • • 
on the supersymmetry parameters \x (left) and tan /3 (right), for —2000 GeV < \i < 2000 GeV 
and 1 .5 < tan /3 < 30, with the others as in Eq. (|3TI) . The distribution of the masses of X\i ■■■ X4 
is symmetric with respect to ji = 0. Nevertheless, m^o and m^o increase from (/1 = 0) to about 
100 (m^o) and 200 GeV (m^o) in the range < 300 GeV, whereas they are almost constant 

for 300 GeV < < 2000GeV. On the contrary, the masses of X3 an d X4 exhibit a minimum 
for /i = 0, about 110 and 230 GeV respectively, and increase monotonically in terms of 
with a behaviour leading to m^o ~ m*o ~ \fl\ for large As for tan/3, a small dependence 

is visible only in the low tan/3 range, i.e. 1.5 < tan/3 < 8, with the masses of %\-> x® an d Z3 
slightly decreasing and the one of £4 mildly increasing. Outside this range, the light neutralino 
masses are roughly independent of tan/3. 

In Fig. [5]we present the dependence of the light (left) and heavy (right) neutralino masses 
on the gaugino mass M\ for M\ < 3.7 TeV. In the light case, the masses exhibit a step-like 
behaviour: m^o and m^o have roughly the same value through all Mi range, growing for small 
Mi and amounting to approximately 200 GeV for Mi > 200 GeV. The mass m^o increases in 
the range 200 GeV < M x < 2.5 TeV and is about m.o ~ 2.54 TeV for M x > 2.5 TeV. The 
mass of X4 is roughly m^o ~ 2M\ for 200 GeV < Mi < 1.2 TeV, then m^o ~ 2.54 TeV, up to 
Mi ~ 2.5 TeV, and ultimately m 9 o ^ Mi for larger Mi. As for the heavy neutralinos, the mass 
of Xs is m v° — 2-54 TeV for Mi < 1.3 TeV, then it increases linearly in the range 1.3 TeV < 

3 A5 

Mi < 1.8 TeV and it is m 9 o — 3.54 TeV for Mi > 1.8 TeV. The mass of the heaviest neutralino 

A5 

X® is constant, namely m^o ~ 3.54 TeV, for Mi < 1.8 TeV, then it grows linearly, reaching the 
value m^o ~ 7 TeV for Mi =3.5 TeV. 

■A6 

Figure |6]presents the masses of X5 an d X® with respect to the Z' mass in the range 1 TeV< 
m z i < 4 TeV (left) and to the M parameter for 100 GeV< M' < 4 TeV (right). The masses 
of x® and %® grow linearly as a function of m z /, whereas they exhibit opposite behaviour with 
respect to M', as m^o increases from 3 to 5.5 TeV and m^o decreases from 3 to 1.5 TeV. The 

A6 X5 

four light-neutralino masses are instead roughly independent of m z i and M', as expected. 
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Figure 4: Dependence of the mass of the four lightest neutralinos on the MSSM parameters /i 
(left) and tan/3 (right). 
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Figure 5: Dependence of the neutralino masses on the MSSM parameter M\. Left: light neu- 
tralinos. Right: heavy neutralinos. 

4.3 Chargino masses 

As discussed before, the chargino sector remains unchanged after the introduction of the extra 
group U(l)'. Therefore, the chargino masses do not depend on the U(l) ; new parameter M' and 
on m z i, but just on the MSSM parameters ji, tan /3 and M\ . Figures|7]and[8]show the dependence 
on such quantities, which are varied individually, with the other parameters fixed as in Eq. (13TT >. 

The dependence on /i, displayed in Fig. [7] (left), is symmetric with respect to /i = 0. In 
particular, m f ± varies significantly, from about 3 to 200 GeV, only for < 300 GeV, whereas 
the heavier chargino mass exhibits a behaviour m 9 ± ~ \a\ and is as large as 2 TeV for lul ~ 

2000 GeV. As for tan/3, Fig. [7] (right), the mass of the heavy chargino %2 increases quite mildly 
from 230 to about 263 GeV, whereas m^± decreases from almost 200 GeV (tan/3 = 1.5) to 
about 154 GeV (tan/3 = 30). 

The variation with respect to M\, presented in Fig. [8l is instead quite different for the two 
charginos. The mass of the lighter one changes very little only for M\ < 200 GeV, whereas for 
larger Mi it is about m^± ~ 200 GeV. The mass of increases almost linearly with M\ and 
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Figure 6: Mass spectra of the heavy neutralinos as functions of the Z' mass (left) and the gaugino 
mass M' (right). 
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Figure 7: Chargino masses with respect to the the MSSM parameters ji (left) and tan /3 (right). 



and is m ± ~ 2M\ ~ M2 for large M\ . 

4.4 Higgs masses 

As pointed out before in the paper, after adding the U(l)' symmetry, one has an extra neutral 
scalar Higgs, named H' , besides the Higgs sector of the MSSM, i.e. the bosons h, H, and 
A. The Z' phenomenology will thus depend on the three Higgs masses and vacuum expectation 
values vi, V2 and V3. In the Representative Point parametrization, the lightest h has a mass 
m/j ~ 90 GeV, H, A and // ± are degenerate and have a mass of about 1 190 GeV, whereas the 
U(l /-inherited H' is about 3 TeV, like the Z' . Therefore, in this scenario the Z' is not capable 
of decaying into final states containing H'. 

Figure [9] presents the variation of the Higgs masses in terms of /1 (left) and tan/3 (right); 
Fig. [TO] shows the dependence on m z i (left) and Ay (right). One can immediately notice that the 
mass of the lightest h is roughly independent of these quantities and it is ~ mz — 90 GeV 
through the whole /i, tan/3, m z t and Ay ranges. Since the supersymmetric light Higgs h should 
roughly play the role of the SM Higgs boson, a value of about 90 GeV for its mass is too 
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Figure 8: Dependence of the chargino masses on the MSSM parameter Mi. 

low, given the current limits from LEP [|36l and Tevatron [37] experiments and the recent LHC 
results [|38~l l39l on the observation of a new Higgs-like particle with a mass about 125 GeV. This 
is due to the fact that the h mass obtained after diagonalizing the neutral Higgs mass matrix is 
just a tree-level result; the possible inclusion of radiative corrections should increase the light 
Higgs mass value in such a way to be consistent with the experimental limits. In fact, the 
Representative Point will be used only to illustrate the features of the particle spectra in the 
MSSM, after one adds the extra U(l)' symmetry group. Any realistic analysis of Z' decays in 
supersymmetry should of course use values of the Higgs masses accounting for higher-order 
corrections and in agreement with the experimental data. 

The heavy MSSM scalar Higgs H is physical, i.e. its squared mass positive definlite, only 
for positive values of /i, therefore in Fig.|9]the Higgs masses are plotted for /i > 0. The mass of 
H increases monotonically from (/i = 0) to 3 TeV (/i ~ 1260 GeV), and then it is m# ~ mz> 
also for larger /i-values. As for the U(l /-inherited H', its mass is about m H > ~ m Z ' for < \l < 
1260 GeV; for larger \x it increases monotonically, up to m H i ~ 3.75 TeV, value reached for }i = 
2000 GeV. In other words, for ji > 1260 GeV, H and H' behave as if they exchanged their roles, 
with increasing m H i and constant m# = m z i. The masses of A and H ± exhibit instead the same 
behaviour and increase monotonically with respect to ji in the whole range. It is also interesting 
to notice that, for < \l < 1260 GeV, one has mn — m H ± ~ m^. As for the dependence on tan/3, 
presented in Fig. [9] (right), the masses of H, A and // ± are almost degenerate and increase from 
about 400 GeV (tan j8 = 1.5) to approximately 1.5 TeV (tan /3 = 30). The mass of H' is instead 
m H i ~ m z i = 3 TeV for any value of tan/3. 

The dependence of the Higgs masses on the Z' mass in the range 1 TeV< m z / < 4 TeV is 
presented in Fig. [101 (left). A and H ± are degenerate and their mass is constantly equal to 1.19 
TeV in the whole explored region. The H mass is m# ~ 1 TeV for m Z ' = 1 TeV, then it slightly 
increases and amounts to m# ~ 1.19 TeV in the range 1.2 TeV< m z / < 4 TeV. Figure ITOl (right) 
shows the Higgs masses as functions of the trilinear coupling Af for 500 GeV < Af < 4 TeV. 
The masses of the charged and pseudoscalar Higgs bosons are degenerate and increase from 
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Figure 9: Dependence of the mass of the Higgs bosons h, H, A, H' and // ± on the MSSM 
quantities /i (left) and tan/3 (right). 
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Figure 10: Higgs mass spectra with respect to the Z' mass (left) and the trilinear coupling A/ 
(right). 

1.1 TeV (A f = 500 GeV) to about 3.4 TeV (A f = 4 TeV). The mass of the scalar neutral H is 
degenerate with the ones of A and H ± for 500 GeV < Af < 3.2 TeV, then it is m# = m z t = 
3 TeV for Af between 3.2 and 4 TeV. The H' mass is constant, i.e. m H i = m z t = 3 TeV for 
500 GeV < Af < 3.2 TeV, then it increases in the same manner as the masses of H ± and A. As 
already observed for the /i dependence, H and H' exchange their roles for Af > 3.2 TeV. 

4.5 Consistency of the MSSM masses with ISAJET 

An experimental search for supersymmetric Z' decays demands the implementation of our 
MSSMAJ(l)' scenario in a Monte Carlo event generator. Therefore, it is essential to verify 
whether our mass spectra are consistent with those provided by the codes typically used to 
compute masses and decay rates in supersymmetry. For this purpose, a widely used program is 
the ISAJET package [40], containing all the MSSM data; the supersymmetric particle masses 
and branching ratios obtained by running ISAJET are then used by programs, such as HER- 
WIG ||2~2l and PYTHIA E3l . which simulate hard scattering, parton showers, hadronization 
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Table 4: Mass values in GeV for neutralinos, charginos and Higgs bosons in our model, based 
on U( 1 )' and the MSSM, and according to the IS AJET code, which implements only the MSSM. 



Model 






It, 


M 


m h 


m H 


m A 


m H ± 


i\ 


m 9 ± 

A? 


U(1)7MSSM 


94.6 


156.6 


212.2 


261.0 


90.7 


1190.0 


1190.0 


1190.0 


155.0 


263.0 


MSSM 


91.3 


152.2 


210.2 


266.7 


114.1 


1190.0 


1197.9 


1200.7 


147.5 


266.8 



and underlying event, for an assigned MSSM configuration. It is thus crucial assessing whether 
such an approach can still be employed even after the inclusion of the extra Z' boson. Squark 
and slepton masses, corrected by the D-term contribution, can be directly given as an input 
to ISAJET. Moreover, the chargino spectrum is unchanged, being the Z' neutral, whereas the 
extra H', being too heavy, is not relevant for the Z' phenomenology. Besides, the masses of 
the MSSM Higgs bosons h, H, A and H ± depend very mildly on the U(l)' parameters. In the 
neutralino sector, the two additional %® an d %® are also too heavy to be phenomenologically 
relevant. However, the neutralino mass matrix, Eq. (|2TT >. depends also on extra new parame- 
ters, such as M', g' and the U(l)' charges Q[ 2 3- Therefore, even the mass of the four light 
neutralinos can potentially feel the effect of the presence of the Z'. 

We quote in Table|4]the eigenvalues of the neutralino mass matrix, Eq. (|2TT >. along with the 
masses yielded by ISAJET, for the parameter configuration corresponding to the Representative 
Point, Eq. (I3TI) . For the sake of completeness, we also present the Higgs and chargino mass 
values obtained in our framework (U(l)' and MSSM), to investigate whether they agree with 
the ISAJET results (only MSSM). 

From Table @] one learns that the masses of the neutralinos agree within 5%; a larger dis- 
crepancy is instead found, about 20%, for the mass of the lightest Higgs, i.e. h; as pointed 
out before, this difference is due to the fact that, unlike ISAJET, our calculation is just a tree- 
level one and does not include radiative corrections. Both Higgs masses are nevertheless much 
smaller than the m z >, fixed to 3 TeV in the Representative Point; therefore, Z' decays into Higgs 
bosons will not be significantly affected by this discrepancy. 

As for the chargino masses, the difference between our analytical calculation and the pre- 
diction of ISAJET is approximately 5% for and 1% for %2- Overall, one can say that 
some differences in the spectra yielded by our computations and ISAJET are visible, but they 
should not have much impact on Z' phenomenology. The implementation of the U(l)' model 
in HERWIG or PYTHIA, along with the employment of a standalone program like ISAJET for 
masses and branching ratios in supersymmetry, may thus provide a useful tool to explore Z' 
phenomenology in an extended MSSM. 

4.6 Z! decays in the Representative Point 

Before concluding this section, we wish to present the branching ratios of the Z' boson into 
both SM and new-physics particles. If BSM decays are competitive with the SM ones, then the 
current limits on the Z' mass will have to be reconsidered. We shall first present the branching 
ratios in the Representative Point parametrization, Eq. (|3T| ), i.e. a Z[ boson with mass 3 TeV, 
and then we will vary the quantities entering in our analysis. 
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4.6.1 Branching ratios in the Representative Point 

In Table[5]we summarize, for the reader's convenience, the masses of the BSM particles for the 
parameters in Eq. (|3TT ). in such a way to figure out the decay channels which are kinematically 
permitted. At this point it is possible to calculate the Z' widths into the kinematically allowed 

Table 5: Masses in GeV of BSM particles in the MSSM/U(1)' scenario, with the parameters set 
as in Eq. (|3TI) . 



mm 


m 2 




"V, 


m l 


m h 


m t 


My, 


2499.4 


2499.7 


2500.7 


1323.1 


3279.0 


2500.4 


3278.1 


3279.1 


















94.6 


156.5 


212.2 


260.9 


2541.4 


3541.4 


154.8 


262 A 


m h 


m A 


m H 


m H i 


m H ± 








90.7 


1190.7 


1190.7 


3000.0 


1193.4 









decay channels. The Z' SM decay channels are the same as the Z boson, i.e. quark or lepton 
pairs, with the addition of the W + W~ mode, which is accessible due to the higher Z' mass. 
However, since the Z' has no direct coupling to W bosons, the Z' — » occurs only via ZZ' 

mixing and therefore one can already foresee small branching ratios. Furthermore, the extended 
MSSM allows Z' decays into squarks, i.e. qiq* (q — u,d and i = 1,2), charged sleptons 
sneutrinos Vi } £Vfi {£ = e, /i, %, i = 1,2), neutralino, chargino, or Higgs (hh, HH, hH, hA, HA, 
H'A, H + H ) pairs, as well as into states with Higgs bosons associated with W/Z, such as Zh, 
ZH and W ± H^. 

We refer to [fl2~l for the analytical form of such widths, at leading order in the U( 1)' coupling 
constant, i.e. &{g' 2 ); in Appendix A the main formulas will be summarized. Summing up all 
partial rates, one can thus obtain the Z' total width and the branching ratios into the allowed 
decay channels. 

In Table [6] we quote the Z' branching ratios in the Representative Point parametrization. 
Since, at the scale of 3 TeV, one does not distinguish the quark or lepton flavour, the quoted 
branching ratios are summed over all possible flavours and uu, dd, £ + £~ and vv denote any 
possible up-, down-type quark, charged-lepton or neutrino pair. Likewise, uu*, dd*, £ + £~ and 
vv* are their supersymmetric counterparts. We present separately the branching ratios into 
all possible different species of charginos and neutralinos, as they yield different decay chains 
and final-state configurations. In Tabled several branching ratios are zero or very small: the 
decays into up-type squarks and sleptons, heavy neutralinos %® and the U(l /-inherited H' are 
kinematically forbidden for a Z' of 3 TeV. The only allowed decay into sfermion pairs is the 
one into down-type squarks ^2^2- Despite being kinematically permitted, the width into up- 
type quarks vanishes, since, as will be clarified in Appendix A, in the Z[ model the vector (v u ) 
and vector-axial (a u ) couplings, contained in the in the interaction Lagrangian of the Z' with up 
quarks, are zero. From Table |6]we learn that, at the Representative Point, the SM decays account 
for roughly the 77% of the total Z' width and the BSM ones for the remaining 23%. As for the 
BSM modes, the rate into down squarks is about 9% of the total rate, the ones into charginos 
and neutralinos 4.2% and 8.4%, respectively. In the gaugino sector, the channels ^2^3 an( ^ 
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Table 6: Branching ratios of the Z ! with the parameters fixed as in Eq. (|3TT ). The branching 
ratios into fermions and sfermions have been summed over all the possible flavours, e.g. uu 
(£ + £~) denotes the sum of the rates into up, charm and top (electron, muon and tau) pairs. 
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X1X2 have me highest branching ratios. The decay into X\X\ has a very small branching 
fraction and is experimentally undetectable if %i is the lightest supersymmetric particle (LSP). 
The final states with Higgs bosons are characterized by very small rates: the branching fractions 
into H + H and HA are about 0.5%, the one into H ± W^ P roughly 0.1% and an even lower rate, 
<^(10~ 7 ), is yielded by the modes hZ, ZH and hA. 

These considerations, obtained in the particular configuration of the Reference Point, Eq. (1311 ), 
can be extended to a more general context. We can then conclude that the Z' BSM branching 
fractions are not negligible and should be taken into account in the evaluation of the mass limits. 

4.6.2 Parameter dependence of the branching ratios 

In this subsection we wish to investigate how the Z' branching fractions into SM and supersym- 
metric particles fare with respect to the U(l)' and MSSM parameters. As in Section 3, the study 
will be carried out at the Representative Point, varying each parameter individually. 

In Fig. [TH the dependence of the branching ratios on the mixing angle 6 is presented for 
SM (left) and BSM (right) decay modes, in the range — 1 < 6 < 0.8; for the SM channels, we 
have also plotted the total branching ratio. The Z' decay rate into quarks exhibits a quite flat 
distribution, amounting to about 40% for central values of and slightly decreasing for large 
1 6 1 . The branching ratio into neutrino pairs is enhanced for 6 at the edges of the explored range, 
being about 25%, and presents a minimum for 6 ~ —0.1. The rate into charged leptons varies 
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between 5 and 15%, with a small enhancement around |0| — 0.8; the branching fraction into 
W + V7 - is below 2% in the whole range. 

As for the BSM channels, described in Fig. [TTI (right), the neutralino, chargino, and Higgs 
modes have a similar behaviour, with a central broad maximum around 6 = and branching 
ratios about 20%, 10% and 3%, respectively The sneutrino modes give a non-negligible contri- 
bution only for > 0.5, reaching about 10%, at the boundary of the investigated region, i.e. 
~ 0.8. The squark-pair channel has a significant rate, about 15%, for negative mixing angles, 
i.e. ~= —1. The rates in the Higgs channels lie between the neutralino and chargino ones 
and exhibit a maximum value, about 10%, for 0=0. 

Figure [[^presents the dependence of the BSM Z' branching ratios on the MSSM parameters 
/l (left) and tan/3 (right). The SM rates are not shown, since their dependence on these param- 
eters is negligible. The decay rate into squarks slightly increases from 9 to 10% in the explored 
jli range; the neutralino branching ratio decreases quite rapidly from about 8% (/i = 0) to zero 
(/l ~ 1500 GeV). The rate into charginos is about 4% for small values of /i, then it smoothly de- 
creases, being negligible for /i > 1500 GeV. The branching fraction into Higgs modes is almost 
4% at /i = and rapidly becomes nearly zero for /i > 300 GeV. As for tan/3, the qq*, X + X~ 
and X°X° modes are roughly independent of it, with rates about 9% (squarks), 8% (neutralinos) 
and 4% (charginos). The decays into states with Higgs bosons account for 4% of the Z' width 
at small tan /3 and are below 1% for tan /3 > 20. 
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Figure 11: Dependence of the Z' decay rates on the U(l)' mixing angle 0. Left: SM modes; 
right: BSM channels. 
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Figure 12: BSM branching ratios with respect to the MSSM parameters /i (left) and tan/3 
(right). 

5 7J decays into final states with leptons 

Leptonic final states are considered as golden channels from the viewpoint of the LHC exper- 
imental searches. To exploit these features, this study will be focused on the decays of the 
Z' boson into supersymmetric particles, leading to final states with leptons and missing en- 
ergy, due to the presence of neutralinos or neutrinos. Final states with two charged leptons and 
missing energy come from primary decays Z' — > £ + £~, presented in Fig. \T3\ with the charged 
sleptons decaying into a lepton and a neutralino. 

Furthermore, primary decays into charginos Z' — > X7X2' followed by %2 ~~ ^ W^Xi an d 
W + — > £ + v (W~ — > £~v), as in Fig. [fj] yield final states with two charged leptons and missing 
energy as well. With respect to the direct production in pp collisions, where the partonic centre- 
of-mass energy is not uniquely determined, the production of charginos in Z' decays has the 
advantage that the Z' mass sets a kinematic constrain on the chargino invariant mass. 

A decay chain, leading to four charged leptons and missing energy, is yielded by Z' decays 
into neutralinos Z' — > X^X®* wmi subsequent %2 — > ^ ± ^ =F an d — > P^Xi* as i n Fig- EES Finally, 
we shall also investigate the decay into sneutrino pairs, such as Z' — > V2V2, followed by V2 — > 
X^y and X2 ~^ Z?' w l tn an intermediate charged slepton (see Fig. [16]). The final state of 
the latest decay chain is made of four charged leptons plus missing energy, due to neutrinos and 
neutralinos. 

In the following, we wish to present a study of Z' decays into leptonic final states for a given 
set of the MSSM and U(l)' parameters. In particular, we shall be interested in understanding 
the behaviour of such rates as a function of the slepton mass, which will be treated as a free 
parameter. In order to increase the rate into sleptons, with respect to the scenario yielded by the 
Representative Point, the squark mass at the Z' scale will be increased to 5 TeV, in such a way 
to suppress Z' decays into hadronic jets. 

In our study we consider the models in Table \T\ and vary the initial slepton mass m~ for 
several fixed values of m z i, with the goal of determining an optimal combination of m~ and m z /, 
enhancing the rates into leptonic final states, i.e. the decay modes containing primary sleptons, 
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charginos or neutralinos. The other parameters are set to the following Reference Point: 

/x = 200 GeV , tan/3 = 20 , A q = A e = A f = 500 GeV , 

m\ = 5 TeV , M x = 150 GeV , M 2 = 300 GeV , M' = 1 TeV. (33) 

Any given parametrization will be taken into account only if the sfermion masses are physical 
after the addition of the D-term. Hereafter, we denote by BR^, BR^ , BR V y and BR^+w the 
branching ratios into quark, charged-lepton, neutrino and W pairs, with BRsm being the total 
SM decay rate. Likewise, BR^*, BR f - +f -- and BRyy* are the rates into squarks, charged sleptons 
and sneutrinos, BR^+^ , BR^o^o, BR H+H ,BR hA , BR H a are the ones into chargino, neutralino, 
charged- and neutral-Higgs pairs, BR W ^ H ± the branching fraction into W^ P H ± . Moreover, for 
convenience, BRzh is the sum of the branching ratios into Zh and ZH and BRbsm the total BSM 
branching ratio. 
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Figure 13: Diagram for the decay of the Z' into a charged- slep ton pair, yielding a final state 
with two charged leptons and missing energy. 
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Figure 14: Final state with two charged leptons and missing energy, through a primary decay of 
the Z' into a chargino pair. 
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Figure 15: Decays of Z' bosons into neutralinos, leading to final states with four charged leptons 
and missing energy. 
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Figure 16: Final state with four charged leptons and missing energy, due to the presence of 
neutrinos and neutralinos, yielded by a primary Z ! — > vv* decay. 



5.1 Reference Point: Model 

An extra U(l)' group with a mixing angle 6 = arccos a/5/8 leads to a new neutral boson la- 
belled as Z^. In Table |7] we list the masses of charged (nig and mg) and neutral (my, and my 2 ) 
sleptons, for various m z > and for the values of m~ which, as will be clarified later, yield a phys- 
ical sfermion spectrum and a maximum and minimum rate into sneutrinos. From Table [7] we 
learn that the decays into pairs of charged sleptons are always kinematically forbidden, whereas 
the decay into V2 pairs is accessible. The effect of the D-term on the mass of V2 is remarkable: 
variations of m~ of few hundreds GeV induce in niy 2 a change of 1 TeV or more, especially 
for large values of the Z' mass. Table [8] summarizes the branching ratios into all allowed SM 
and BSM channels, for the same and m~ values as in Table [71 whereas Fig. [F7]presents the 
branching ratio Z^ — > V2V2 as a function of ml and for 1 TeV < m z i < 4 TeV. The branching 
fraction into sneutrinos can be as large as about 11% for any value of m z r, for larger m~ the 
sneutrino rate decreases, as displayed in Fig. [FT] Furthermore, Table [8] shows that, within the 
scenario identified by the Reference Point, even the decays into charginos and neutralinos are 
accessible, with branching ratios about 5-6% (charginos) and 10-12% (neutralinos). Decays 
into VK + VK~ pairs or Higgs bosons associated with Z's are also permitted, with rates about 3%. 
The decrease of the sneutrino rate for large ml results in an enhancement of the SM branching 
ratios into qq and neutrino pairs. As a whole, summing up the contributions from sneutrinos, 
charginos and neutralinos, the branching ratio into BSM particles runs from 24 to 33%, thus 
displaying the relevance of those decays in any analysis on Z' production in a supersymmetric 
scenario. 
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Figure 17: Branching ratio of the boson into sneutrino pairs V2V2, as a function of the slepton 
mass ml, for several values of m^'- 



Table 7: Slepton masses at the Reference Point for the Z^ model, varying m z / and m~. 2 
and my l 2 are the charged- slepton and sneutrino mass eigenvalues, as discussed in the text. All 
masses are given in GeV. 
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Table 8: Branching ratios of the boson into SM and BSM channels, varying m z / and nv~, 
given in TeV, along the lines described in the text. BRsm and BRbsm denote the total branching 
fractions, respectively. 
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5.2 Reference Point: Zy 

An extra group U(l)' with a mixing angle 0=0 leads to a neutral vector boson labelled as Z' 
(Tabled]). In Tabled we quote the slepton masses for a few values of m z i and nf~: as before, 
the results are presented for the two values of m~ which are found to enhance and minimize 
the slepton rate. For any mass value, the D-term enhances by few hundreds GeV the masses 
of l\ and Vi and strongly decreases mg 2 and wzj , especially for small m~ and large m z /. In 
Table [TOl we present the branching ratios into all channels, for the same values of m z > and as 
in Tabled Unlike the Z^ case, supersymmetric decays into charged- slepton pairs are allowed 
for = 0, with a branching ratio, about 2%, roughly equal to the sneutrino rate. Furthermore, 
even the decays into gauginos are relevant, with rates into X X an d X°X° about 10 and 20%, 
respectively. The decays into boson pairs, i.e. Zh and W + V7~, are also non-negligible and 
account for about 3% of the total Z^ width. 

As a whole, the Z'^ modelling above depicted yields branching ratios of the order of 35- 
40% into BSM particle, and therefore it looks like being a promising scenario to investigate 
Z' production within the MSSM. Figure [TBI finally displays the branching ratios into sneutrinos 
and charged sleptons as a function of m~ and for several values of mz 1 - 

5.3 Reference Point: 

In this subsection we investigate the phenomenology of the Z' N boson, i.e. a U(l)' gauge group 
with a mixing angle = arctan \/T5 — n/2 (Tabled)), along the lines of the previous sections. As 
discussed above, the Z' N model is interesting since it corresponds to the Z' model, but with the 
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Table 9: Slepton masses in GeV at the Reference Point for the model Z'y and a few values of 
m z i and m~. 







m h 


m h 


ffiv, 


My, 


1000 


400 


535.2 


194.2 


529.2 


189.2 


1000 


700 


785.1 


606.4 


781.0 


604.8 


1500 


600 


801.7 


285.4 


797.7 


282.0 


1500 


1000 


1132.6 


849.4 


112.7 


848.3 


2000 


800 


1068.4 


377.8 


1065.4 


375.2 


2000 


1300 


1480.3 


1092.1 


1478.2 


1091.2 


2500 


1000 


1335.2 


470.6 


1333.8 


468.6 


2500 


1600 


1828.3 


1334.7 


1826.6 


1334.0 


3000 


1100 


1528.5 


296.2 


1526.4 


292.9 


3000 


1900 


2176.3 


1577.2 


2174.9 


1576.6 


3500 


1300 


1795.2 


401.8 


1793.4 


399.4 


3500 


2200 


2524.4 


1819.7 


2523.2 


1819.2 


4000 


1500 


2061.9 


502.7 


2060.4 


500.8 


4000 


2500 


2872.5 


2062.2 


2871.4 


2061.7 


4500 


1600 


2256.7 


177.4 


2255.3 


171.9 


4500 


2800 


3220.7 


2304.7 


3219.7 


2304.2 


5000 


1800 


2523.2 


343.1 


2521.9 


340.3 


5000 


3100 


3568.8 


2547.1 


3567.9 


2546.7 




m,,=5.0TcV 

— m z ,=4.5TeV 
.... m,=4.0TeVT 
.... m z ,=3.5TeV : 

— m z .=3.0TeV 
.... m z =2.5TeV 

— m z ,=2.0TeV 

— m z ,=l.STeV 
... m..=1.0TeV 



500 



1000 1500 2000 2500 3000 3500 4000 

m»[GeV] 



•| °- 02 
1 0.018 

DC 

J 0.016 
1 0.014 
« 0.012 
0.01 
0.008 
0.006 
0.004 
0.002 





— m z .=4.5TeV 
.... m z .=4.0TeV^ 
... m z .=3.STeV 

m z .=3.0TeV 

..... m z .=2.5TeV 

— m z .=2.0TeV 

— m z .=l.STeV 
. . . m,.=1.0 TeV 



500 1000 1500 2000 2500 3000 3500 4000 

m»[GeV] 



Figure 18: Dependence of the branching ratio into charged sleptons (left) and sneutrinos 
(right) as a function of rr&, for several values of m z i. 
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Table 10: Branching ratios of the boson into SM and BSM channels, varying m z i and m~. 
The masses are expressed in TeV. 
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unconventional assignment of the SO(10) representations. Referring to the notation in Eq. ©, 
in the unconventional Eg model the fields H and D c are in the representation 16 and L and d c in 
the 10 of SO(10). 

Table ITTI presents the slepton masses varying m z / and for the values of m~ which minimize 

and maximize the slepton rate. The D-term addition to m~ increases the mass of l\ and Vi 

and decreases the mass of £2; its impact on V2 is negligible and one can assume my 2 — m~. 

Both decays into t\i 2 and V2V2 316 kinematically allowed, whereas l\ and Vi are too heavy to 
contribute to the Z' width. 

Table [T2l quotes the branching ratios for the Z^, computed for the same values of and 
m~ as in Table QT] Although — > V2V| is kinematically allowed, the coupling of the Z^ to 

sneutrinos is zero for = arctan vT5 — Tt/2, since, as will be discussed in Appendix A, the rate 
into right-handed sfermions vanishes for equal vector and vector-axial coupling, i.e. Vy = ay. 
therefore, this decay mode can be discarded. As for the other supersymmetric channels, the 
rates into charginos and neutralinos are quite significant and amount to about 9% and 28%, 
respectively. The decays into V7 + V7~ and Zh states account for approximately 1-2%, whereas 
the branching ratio into charged- slepton pairs is about 1%, even in the most favourable case. As 
a whole, the rates into BSM final states run from 18 to about 35% and therefore are a relevant 
contribution to the total Z' cross section. Figure [19] finally presents the variation of the charged- 
slepton branching ratio as a function of m~, for a few values of m Z '. 

Table 1 1 : Slepton masses in the Z^ model, varying m Z ' and m?, as discussed in the text. All 
masses are given in GeV. 



m z > 




m h 


m h 






1000 


400 


601.1 


249.7 


595.8 


400.0 


1000 


600 


749.2 


512.2 


745.0 


600.0 


1500 


500 


837.4 


165.4 


833.6 


500.0 


1500 


900 


1123.1 


766.4 


1120.2 


900.0 


2000 


700 


1136.4 


303.9 


1133.6 


700.0 


2000 


1200 


1497.1 


1021.0 


1495.0 


1200.0 


2500 


800 


1375.8 


131.8 


1372.9 


800.0 


2500 


1500 


1871.2 


1275.7 


1869.5 


1500.0 


3000 


1000 


1673.7 


319.9 


1671.8 


1000.0 


3000 


1800 


2245.3 


1530.4 


2243.9 


1800.0 


3500 


1200 


1972.6 


466.2 


1971.0 


1200.0 


3500 


2100 


2619.4 


1785.3 


2618.2 


2100.0 


4000 


1300 


2211.6 


303.9 


2210.2 


1300.0 


4000 


2400 


2993.6 


2040.2 


2992.5 


2400.0 


4500 


1500 


2510.2 


476.8 


2509.0 


1500.0 


4500 


2700 


3367.7 


2295.1 


3366.7 


2700.0 


5000 


1600 


2749.8 


249.7 


2748.6 


1600.0 


5000 


3100 


3822.5 


2666.9 


3821.6 


3100.0 



32 



Table 12: Branching ratios of the boson in SM and BSM channels, varying m z i and m*. 
Slepton and Z' masses are quoted in TeV. 



mz 1 


m l 


RR - 


BR(5+/;- 


RR - 


T3T? 




RR _l 


RR n n 

XX 


RR - - 


BRsm 


CK BSM 


1.0 


0.4 


49.51 


11.98 


9.59 


1.71 


1.68 


8.71 


15.78 


1.04 


72.79 


27.21 


1.0 


0.6 


50.03 


12.11 


9.69 


1.73 


1.69 


8.80 


15.94 


0.00 


73.56 


26.44 


1.5 


0.5 


47.99 


11.51 


9.21 


1.57 


1.57 


9.26 


17.76 


1.12 


70.28 


29.72 


1.5 


0.9 


48.53 


11.64 


9.31 


1.59 


1.59 


9.36 


17.96 


0.00 


71.08 


28.92 


2.0 


0.7 


47.50 


11.36 


9.08 


1.53 


1.54 


9.44 


18.46 


1.08 


69.47 


30.53 


2.0 


1.2 


48.02 


11 AO 

11.48 


9.18 


1 C A 

1.54 


1 cc 

1.55 


9.55 


1 O S f 

18.66 


0.00 


70.22 


or* to 

29.78 


2.5 


0.8 


47.16 


11.26 


9.01 


1.50 


1.52 


9.50 


18.73 


1.12 


68.92 


31.08 


2.5 


1.5 


47.69 


11.38 


9.11 


1.52 


1.53 


9.61 


18.94 


0.00 


69.70 


30.30 


3.0 


l.O 


46.43 


11.30 


8.86 


1.47 


1.49 


9.43 


18.66 


1.08 


67.83 


32.17 


3.0 


1.8 


46.94 


11.20 


8.96 


1.49 


1.50 


9.53 


18.86 


0.00 


68.58 


31.42 


3.5 


1.2 


45.85 


10.93 


8.74 


1.45 


1.47 


9.35 


18.56 


1.05 


66.98 


33.02 


3.5 


2.1 


46.34 


11.05 


8.84 


1.46 


1.48 


9.45 


18.76 


0.00 


67.68 


32.32 


4.0 


1.3 


45.42 


10.83 


8.66 


1.43 


1.45 


9.29 


18.47 


1.07 


66.34 


33.66 


4.0 


2.4 


45.91 


10.94 


8.75 


1.45 


1.47 


9.39 


18.67 


0.00 


67.06 


32.94 


4.5 


1.5 


45.13 


10.75 


8.60 


1.42 


1.44 


9.24 


18.41 


1.05 


65.90 


34.10 


4.5 


2.7 


45.60 


10.87 


8.70 


1.44 


1.46 


9.34 


18.60 


0.00 


66.61 


33.39 


5.0 


1.6 


44.90 


10.70 


8.56 


l.4l 


1.43 


9.21 


18.35 


1.06 


65.56 


34.44 


5.0 


3.1 


45.38 


10.81 


8.65 


1.43 


1.45 


9.31 


18.55 


0.00 


66.27 


33.73 
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Figure 19: Slepton branching ratios of the boson as a function of m~. 
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5.4 Reference Point: Zj 

The U(l /-based model leading to a Z[, i.e. a mixing angle = arccos a/5/8 — Tt/2, has been 
extensively discussed, as it corresponds to the Representative Point. It exhibits the property that 
the initial slepton mass m~ can be as low as a few GeV, still preserving a physical scenario for 
the sfermion masses. In the following, we shall assume a lower limit of m~ = 200 GeV and 
present results also for 1 TeV, in order to give an estimate of the dependence on m~. 

In Table [T3l the charged- and neutral-slepton masses are listed for a few values of m Z ' and 
m~. We already noticed, when discussing the Representative Point and Fig. [3l that the D-term 
correction to the slepton mass is quite important for l\, V\ and V2, especially for small values 
of m~: this behaviour is confirmed by Table [T3l The D-term turns out to be positive and quite 

large and the only kinematically permitted decay into sfermions is Z[ — > £2^2- However, as in 
the case, the vector and vector-axial coupling are equal, i.e. vg = ag, thus preventing this 
decay mode for the reasons which will be clarified in Appendix A. The conclusion is that in 
the Reference Point scenario, the Z[ boson can decay into neither charged nor neutral sleptons. 
Therefore, the dependence of the branching ratios on m~ is not interesting and Table [141 just 
reports the decay rates for fixed m - = 1 TeV. The total BSM branching ratio lies between 12 
and 17% and is mostly due to decays into chargino (~ 4%) and neutralino (~ 8-9%) pairs. 
Decays involving supersymmetric Higgs bosons, such as H + H , W ± H^ P and HA final states, 
are possible, but with a total branching ratio which is negligible for small Z' masses and at most 
3% for m z > > 4 TeV. As for the decay into SM quarks, it was already pointed out in Table [6] 
that the rate into uu pairs is zero since the couplings v u and a u (see also Appendix A) vanish 
for 6 = arccos a/5/8 — n/2. Therefore, in Table [TU BR^ only accounts for decays into down 
quarks. 

5.5 Reference Point: Z' s 

The Zg boson corresponds to a a mixing angle = arctan(vT5/9) — n/2. As in the Z[ model, 
one can set a small value of the initial slepton mass, such as m~ = 200 GeV, and still have 
a meaningful supersymmetric spectrum. The results on slepton masses and branching ratios 
are summarized in Tables [T31 and [T6] Since the Z' s decay rates are roughly independent of the 
slepton mass, in Table [TBI the branching ratios are quoted only for m~ = 200 GeV. From Table [TBI 
we learn that the D-term contribution to slepton masses is positive and that Z' s — > is the only 
decay kinematically allowed, at least for relatively small values of m~. However, as displayed 
in Table [T5l the branching ratio into such charged sleptons is very small, about 0.1%, even for 
low m? values. As for the other BSM decay modes, the most relevant ones are into chargino 
(about 3%) and neutralino (about 6-7%) pairs, the others being quite negligible. It is interesting, 
however, noticing that for m Z ' = 5 TeV the branching ratio into squark pairs starts to play a role, 
amounting to roughly 8%. In fact, although we set a high value like m~ = 5 TeV, for relatively 

large Z ! masses, i.e. m z / > 3.8 TeV, the D-term for J2-type squarks starts to be negative, in 
such a way that d,2d\ final states are kinematically permitted. As a whole, one can say that, 
at the Reference Point, for m z / < 5 TeV the BSM decay rate is about 10-12%, but it becomes 
much higher for larger Z' masses, even above 20%, due to the opening of the decay into squark 
pairs. However, since the experimental signature of squark production is given by jets in the 
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Table 13: Slepton masses in the U(l)' scenario corresponding to a Z[ boson, for a few values of 
m Z ' and rrR. All masses are expressed in GeV. 







m~ lx mi 2 m Vl m Vl 


1000 
1000 


200 
1000 


736.3 204.7 732.0 734.8 
1226.6 1001.0 1223.0 1224.7 


1500 
1500 


200 
1000 


1080.4 204.7 1077.4 1079.3 

1458.5 1001.0 1456.3 1457.7 


2000 
2000 


200 
1000 


1429.1 204.7 1426.8 1428.3 
1732.7 1001.0 1730.8 1732.0 


2500 
2500 


200 
3000 


1779.7 204.7 1777.9 1779.0 
3482.4 3000.3 3481.5 3482.1 


3000 
3000 


200 
3000 


2131.5 204.7 2129.7 2130.7 
3674.5 3000.3 3673.7 3674.2 


3500 
3500 


200 
3000 


2483.4 204.7 2482.1 2482.9 
3889.4 3000.3 3888.5 3889.1 


4000 
4000 


200 
3000 


2836.9 204.7 2834.8 2835.5 
4123.4 3000.3 4122.6 4123.1 


4500 
4500 


200 
3000 


3188.6 204.7 3187.6 3188.3 
4373.5 3000.3 4372.7 4373.2 


5000 
5000 


200 
3000 


3541.5 204.7 3540.6 3541.2 
4637.0 3000.3 4636.4 4636.8 



Table 14: Branching ratios of the Z[ into SM and BSM particles for m~ = 1 TeV and varying 
m z i. The Z' mass is expressed in TeV. 



m Z ' 




BR^ 


BR^+£- 


BR V y 


BR#+# 


BR W t#± 


BR//A 


BR rr 


BR^o^o 


BRsm 


BRbsm 


1.0 


1.0 


44.06 


14.69 


29.37 


0.00 


£?(10~ 3 ) 


^(io- 4 ) 


4.31 


7.58 


88.11 


11.89 


1.5 


1.0 


43.39 


14.46 


28.93 


0.00 


^(10~ 4 


^(io- 4 ) 


4.56 


8.65 


86.78 


13.22 


2.0 


1.0 


43.16 


14.38 


28.77 


0.00 


^(io- 4 ) 


^(io- 3 ) 


4.65 


9.03 


86.31 


13.69 


2.5 


1.0 


42.99 


14.33 


28.66 


0.06 


£?(10- 3 ) 


0.07 


4.68 


9.19 


85.98 


14.02 


3.0 


1.0 


42.53 


14.18 


28.36 


0.53 


^(10~ 3 ) 


0.53 


4.66 


9.20 


85.07 


14.93 


3.5 


1.0 


42.16 


14.05 


28.11 


0.91 




0.92 


4.64 


9.19 


84.33 


15.67 


4.0 


1.0 


41.90 


13.96 


27.93 


1.20 


^(io- 3 ) 


1.21 


4.62 


9.17 


83.79 


16.21 


4.5 


1.0 


41.70 


13.90 


27.80 


1.40 


£?(10- 3 ) 


1.41 


4.61 


9.16 


83.40 


16.60 


5.0 


1.0 


41.56 


13.85 


27.71 


1.56 


0.01 


1.57 


4.60 


9.15 


83.12 


16.88 
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final state, it is quite difficult separating them from the QCD backgrounds. This scenario seems 
therefore not very promising for a possible discovery of supersymmetry via Z 1 decays. 



Table 15: Slepton masses at the Reference Point with a Z' s gauge boson and for few values of 
mz 1 and m~, given in GeV, as debated in the text. 





4 


m~ lx m- ti my! fflji. 


1 f\f\f\ 

1000 
1000 


200 
1000 


917.9 376.8 914.4 1020.0 
1342.6 1049.7 1340.2 1414.3 


1 C f\f\ 

1500 
1500 


200 
1000 


1357.4 516.7 1355.0 1513.4 
1674.1 1107.7 1672.2 1802.9 


2000 
2000 


200 
1000 


1 O f\f\ 'H ZTZT/IO i nnn (\ 1A1 A A 

1800.7 664.8 1798.9 2010.0 
2050.0 1184.0 2048.4 2236.1 


2500 
2500 


200 
3000 


2245.5 816.7 2244.1 2508.0 
3742.0 3102.7 3741.1 3905.2 


3000 
3000 


200 
3000 


2691.2 970.5 2690.0 3006.7 
4025.2 3146.7 4024.4 4242.7 


3500 
3500 


200 
3000 


3137.3 1125.6 3136.3 3505.7 
4336.2 3198.0 4335.4 4609.8 


4000 
4000 


200 
3000 


3583.6 1281.4 3582.7 4005.0 
4669.3 3256.1 4668.6 5000.0 


4500 
4500 


200 
3000 


4030.2 1437.7 4029.4 4504.5 
5020.2 3320.7 5019.6 5408.4 


5000 
5000 


200 
3000 


4476.9 1594.3 4476.2 5004.0 
5385.4 3391.4 5384.8 5831.0 



5.6 Reference Point: 

The U(l)' group corresponding to a mixing angle 9 = —n/2 and a boson Z' x does not lead to 
a meaningful sfermion scenario in the explored range of parameters, as the sfermion masses 
are unphysical after the addition of the D-term. This feature of the Z' x model, already observed 
in the Representative Point parametrization (see Section 4.1 and the mj 2 spectrum in Fig. [D, 
holds even for a higher initial squark mass, such as m~ = 5 TeV, as in the Reference Point. It 
is nevertheless worthwhile presenting in Table [P71 the Standard Model branching ratios, along 
with those into Higgs and vector bosons in a generic Two Higgs Doublet Model. For any mz' 
the rates into quark and neutrino pairs are the dominant ones, being about 40-45%, whereas the 
branching ratio into lepton states is approximately 12% and the other modes (W + W~, Zh, HA 
and H + H ) account for the remaining 1-3%. 

5.7 Reference Point: Z5 SM 

A widely used model in the analyses of the experimental data is the Sequential Standard Model 
(SSM): in this framework the Z' coupling to SM and MSSM particles is the same as the Z 
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Table 16: Branching ratios of the Z' s with the MSSM parameters at the Reference Point and for 
a few values of m z t, expressed in TeV. The initial slepton mass is fixed to 0.2 TeV, since the 
decay rates are independent of m~. 



m z > 


m\ 


BR^ 


BR(+(i- 


BR V y 


BR^/+^- 


BRz/i 


BR rr 




BR|+£ 


BR^* 


BRsm 


BRbsm 


1.0 


0.2 


42.29 


13.70 


34.57 


0.15 


0.14 


3.33 


5.75 


0.07 


0.00 


90.71 


9.29 


1.5 


0.2 


41.84 


13.54 


34.16 


0.15 


0.14 


3.51 


6.59 


0.07 


0.00 


89.68 


10.32 


2.0 


0.2 


41.67 


13.48 


34.02 


0.14 


0.14 


3.57 


6.90 


0.08 


0.00 


89.32 


10.68 


2.5 


0.2 


41.56 


13.44 


33.91 


0.14 


0.14 


3.59 


7.03 


0.08 


0.00 


89.06 


10.94 


3.0 


0.2 


41.25 


13.34 


33.66 


0.14 


0.14 


3.58 


7.06 


0.08 


0.00 


88.39 


11.61 


3.5 


0.2 


40.99 


13.26 


33.45 


0.14 


0.14 


3.57 


7.07 


0.08 


0.00 


87.84 


12.16 


4.0 


0.2 


40.81 


13.20 


33.30 


0.14 


0.14 


3.56 


7.07 


0.08 


0.00 


87.44 


12.56 


4.5 


0.2 


40.67 


13.15 


33.19 


0.14 


0.14 


3.56 


7.07 


0.08 


0.00 


87.15 


12.85 


5.0 


0.2 


37.34 


12.07 


30.46 


0.13 


0.13 


3.27 


6.50 


0.07 


7.97 


80.00 


20.00 



Table 17: Branching ratios of the Z' x boson as a function of the Z' mass, given in TeV. The rates 
into sfermion pairs are not presented, since the sfermion mass spectrum is unphysical for the 
Z' model in the Reference Point scenario. 



mz 1 


BR^ 


BR^+£ 


BRyy 


BRiy+^- 






BRm 


BRsm 


BRbsm 


1.0 


44.35 


12.44 


42.29 


0.90 


0.00 


0.02 




99.98 


0.02 


2.0 


44.32 


12.34 


41.96 


0.84 


0.00 


0.28 


0.26 


99.46 


0.54 


3.0 


44.03 


12.24 


41.63 


0.82 


0.24 


0.53 


0.52 


98.71 


1.29 


4.0 


43.84 


12.18 


41.43 


0.82 


0.46 


0.64 


0.63 


98.27 


1.73 


5.0 


43.74 


12.15 


41.33 


0.81 


0.58 


0.70 


0.69 


98.03 


1.97 
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boson. The SSM is considered as a benchmark, since the production cross section is only 
function of the Z ! mass and there is no dependence on the mixing angle 6 and possible new 
physics parameters, such as the MSSM ones. 

As for the supersymmetric sector, the sfermion masses get the D-term contribution asso- 
ciated with the hyperfine splitting, Eq. (|25l) . but not the one due to further extensions of the 
MSSM, namely Eq. (f26l) . proportional to g' 2 in the case of U(l)'. Moreover, the Z' SSM cou- 
pling to sfermions is simply given by gssM = gl/ (2cos Qw), as in the SM. Since the hyperfine- 
splitting D-term is quite small, the sfermion spectrum is physical even for low values of m~. 
Table [TBI reports the sfermion masses obtained at the Reference Point, Eq. d33l) . for a few val- 
ues of m z > and varying ml from 100 GeV to m z //2, the highest value kinematically allowed. 
For m~ = 100 GeV, because of the D-term, decreases by about 25%, and slightly 
increase and m.y 2 is roughly unchanged. For large values of m~, the D-term is negligible and all 
slepton masses are approximately equal to m~. 

Tables [191 and 1201 present, respectively, the SM and BSM branching ratios of the Z' SSM at 
the Reference Point, for the values of Z' and slepton masses listed in Table [T8j The decays 
into BSM particles exhibit rates, about 60-65%, which can be even higher than the SM ones, 
accounting for the remaining 35-40%. In fact, this turns out to be mostly due to the decays into 
neutralinos, accounting for more than 30%, and into charginos, about 16-18%. The branching 
fractions into sleptons are quite small: the one into sneutrinos is less than 4% and the one into 
charges sleptons about 1-2%. The V7 + V7~ mode contributes with a rate about 4-5%, the H + H~ 
one is relevant only for m z i > 2.5 TeV, with a branching ratio which can reach 3%, the Zh 
and hA channels are accessible for m z / > 1.5 TeV, with decay fractions between 1 and 4%. The 
variation of the sneutrino and charged- slepton branching ratios as a function of the slepton mass 
at the Z' scale is displayed in Fig.[20]for 1 TeV< m z / <4 TeV. 




m°[GeV] m»[GeV] 

Figure 20: Branching ratios of the Z' SSM as a function of m~ for several values of the Z' mass. 
Left: branching fraction into charged sleptons. Right: branching fraction into sneutrinos. 
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Table 18: Slepton masses in the Z' SSM model varying m Z ' and m~. All masses are quoted in GeV. 



m z > 


,n«, 


m h 


m h 


my, 


m 2 


1000 


100 


110.6 


109.1 


76.6 


100.0 


1000 


500 


502.2 


501.9 


495.8 


500 


1500 


100 


110.6 


109.1 


76.6 


100.0 


1500 


750 


751.5 


751.3 


747.2 


750.0 


2000 


100 


110.6 


109.1 


76.6 


100.0 


2000 


1000 


1001.1 


1000.9 


997.9 


1000.0 


2500 


100 


110.6 


109.1 


76.6 


100.0 


2500 


1250 


1250.9 


1250.8 


1248.3 


1250.0 


3000 


100 


110.6 


109.1 


76.6 


100.0 


3000 


1500 


1001.1 


1000.9 


997.9 


1000.0 


3500 


100 


110.6 


109.1 


76.6 


100.0 


3500 


1750 


1750.6 


1750.6 


1748.8 


1750.0 


4000 


100 


110.6 


109.1 


76.6 


100.0 


4000 


2000 


2000.6 


2000.5 


1999.0 


2000.0 


4500 


100 


110.6 


109.1 


76.6 


100.0 


4500 


2250 


2250.5 


2250.4 


2249.1 


2250.0 


5000 


100 


110.6 


109.1 


76.6 


100.0 


5000 


2500 


2500.4 


2500.4 


2499.2 


2500.0 
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Table 19: Branching ratios into SM particles of the Z' SSM , varying m z i and m®, as debated in the 
text. Slepton and Z' masses are expressed in TeV. 



m z > 




BR^ 




BR V v 


BR w+lv - 


BRsm 


1.0 


0.10 


29.61 


3.87 


7.69 


5.56 


46.73 


1.0 


0.50 


31.38 


4.10 


8.15 


5.90 


49.53 


1.5 


0.10 


27.38 


3.53 


7.02 


4.86 


42.79 


1.5 


0.75 


28.89 


3.73 


7.41 


5.13 


45.15 


2.0 


0.10 


26.21 


3.36 


6.69 


4.56 


40.83 


2.0 


1.00 


27.59 


3.54 


7.04 


4.80 


42.98 


2.5 


0.10 


25.35 


3.25 


6.46 


4.37 


39.42 


2.5 


1.25 


26.64 


3.41 


6.79 


4.59 


41.42 


3.0 


0.10 


24.78 


3.17 


6.31 


4.25 


38.51 


3.0 


1.50 


26.01 


1.66 


6.62 


4.46 


40.42 


3.5 


0.10 


24.42 


3.12 


6.21 


4.17 


37.92 


3.5 


1.75 


25.61 


1.40 


6.51 


4.37 


39.78 


4.0 


0.10 


24.18 


3.09 


6.15 


4.12 


37.54 


4.0 


2.00 


25.35 


1.21 


6.44 


4.32 


39.35 


4.5 


0.10 


24.01 


3.07 


6.10 


4.09 


37.27 


4.5 


2.25 


25.16 


1.07 


6.39 


4.28 


39.06 


5.0 


0.10 


23.89 


3.05 


6.07 


4.06 


37.07 


5.0 


2.50 


25.03 


0.96 


6.36 


4.25 


38.84 



40 



Table 20: Branching ratios into BSM particles of the Zg SM for a few values of m z i and ra~, 
expressed in TeV. 









BRza 










BRvv* 


BRbsm 


1.0 


0.10 


0.00 


~ 10~ b 


0.00 


18.31 


29.30 


1.89 


3.77 


53.27 


1.0 


0.50 


0.00 


~ io~ 6 


0.00 


19.41 


31.06 


0.00 


0.00 


50.47 


1.5 


0.10 


0.00 


0.87 


0.76 


17.84 


32.52 


1.75 


3.48 


57.21 


1.5 


0.75 


0.00 


0.92 


0.80 


18.82 


34.31 


0.00 


0.00 


54.55 


2.0 


0.10 


0.00 


1.93 


1.85 


17.37 


33.01 


1.67 


3.33 


59.17 


2.0 


1.00 


0.00 


2.04 


1.95 


18.28 


34.75 


0.00 


0.00 


57.02 


2.5 


0.10 


0.91 


2.59 


2.53 


16.93 


32.78 


1.62 


3.22 


60.58 


2.5 


1.25 


0.95 


2.72 


2.66 


17.79 


34.45 


0.00 


0.00 


58.57 


3.0 


0.10 


1.72 


2.98 


2.94 


16.62 


32.51 


1.58 


3.15 


61.49 


3.0 


1.50 


1.81 


3.13 


3.08 


17.44 


34.12 


0.00 


0.00 


59.58 


3.5 


0.10 


2.27 


3.23 


3.20 


16.42 


32.30 


1.56 


3.10 


62.08 


3.5 


1.75 


2.38 


3.38 


3.35 


17.22 


33.88 


0.00 


0.00 


60.22 


4.0 


0.10 


2.65 


3.39 


3.37 


16.28 


32.16 


1.54 


3.07 


62.46 


4.0 


2.00 


2.78 


3.56 


3.53 


17.07 


33.71 


0.00 


0.00 


60.65 


4.5 


0.10 


2.91 


3.51 


3.49 


16.19 


32.06 


1.53 


3.05 


62.73 


4.5 


2.25 


3.05 


3.67 


3.65 


16.96 


33.59 


0.00 


0.00 


60.94 


5.0 


0.10 


3.11 


3.59 


3.57 


16.12 


31.98 


1.52 


3.03 


62.93 


5.0 


2.50 


3.26 


3.76 


3.74 


16.89 


33.51 


0.00 


0.00 


61.16 
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6 Cross sections and event rates at the LHC 



In this section we present the total cross section for Z' production at the LHC according to the 
models discussed throughout this paper, i.e. Tabled] as well as the Sequential Standard Model. 
We consider pp collisions at three centre-of-mass energies: 7 TeV (the 201 1 LHC run), 8 TeV 
(the 2012 run) and 14 TeV, the ultimate project energy. For each energy we shall calculate 
the cross section and estimate the expected number of events with a Z' boson decaying into 
supersymmetric particles, for a few values of integrated luminosity. 

6.1 Leading order Z! production cross section 

The cross sections are computed at leading order (LO), employing the LO parton distribution 
functions CTEQ6L |[4TT| and setting the factorization scale equal to the Z' mass. Using a dif- 
ferent LO PDF has a negligible impact on the results. The cross section for Drell-Yan like 
processes has been computed up to next-to-next-to leading order (NNLO) in QCD and, in prin- 
ciple, the calculations carried out in Refs. 114211431 can be easily extended to Z' production 
processes. However, since all Z' partial widths and branching ratios have been evaluated at LO, 
for the sake of consistency, we decided to stick to the lowest-level approximation. 

The parton-level process is analogous to Z production, i.e. it is the purely SM quark- 
antiquark annihilation qq — >• Z ' . Since the coupling of the Z' to the quarks depends on the 
specific U(l)' scenario, the production rate is a function of the mixing angle 6 and of the Z' 
mass, but is independent of the MSSM parameters. In the Sequential Standard Model, the cross 
section just depends on the mass of the Z ; SSM . Figures I2D - I231 present the total cross section for 
the different models investigated throughout this work, as a function of m z /, at the energies of 
7 TeV (Fig. |2~TT) . 8 TeV (Fig. l22l) and 14 TeV (Fig. |23~T) . For each centre-of-mass energy, we 
present the results on linear (left) and logarithmic (right) scales. Tables [2T1 l22l and 1231 quote the 
numerical values of the LO Z' production cross section, varying m z i from 1 to 5 TeV, with steps 
of 500 GeV, in U(l)' models and in the Sequential Standard Model. 

The highest production cross section is given by the SSM, whereas the Z'y model yields 
the lowest rate; the predictions of the other models lie between these results and are almost 
indistinguishable for large m z i. Moreover, the rates decrease by several orders of magnitude 
once m z i increases. In detail, at y/s = 7 TeV, the SSM cross section runs from 1.6 pb (m z / = 
1 TeV) to <^(10^ 8 ) pb (m z t = 5 TeV). The production rate for the U(l /-based Z' varies from 
^(lO -1 ) to ^?(10~ 9 ) pb in the same m Z / range, with very little differences among the models. 
At the centre-of-mass energy of 8 TeV, the variation is between 2.3 pb (Z ; SSM at m z i = 1 TeV) 
and G{ 10~ 9 ) pb (all other models at m z > = 5 TeV). At y/s = 14 TeV, for a Z' mass of 1 TeV the 
cross section varies from about 8 pb (Z' SSM ) to 1.8 pb (ZL); for m z i = 5 TeV, all models yield a 
rate around ^(10~ 4 ) pb. 

6.2 Event rates with sparticle production in Z! decays at the LHC 

In the following, we wish to investigate the domain where possible Z' decays into supersym- 
metric particles could be detectable. For this purpose, we consider two scenarios: ^=8 TeV, 
with an integrated luminosity, / df=20 fb _1 , as expected in the 2012 LHC data taking, and, in 
future perspective, \fs=\A TeV with / Jzf d?=100 fb _1 . In the narrow-width approximation, the 
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Figure 21: Cross section of Z' production in pp collisions at 7 TeV. Left: linear scale. Right: 
logarithmic scale. 




Figure 22: Cross section for Z' production in pp collisions at 8 TeV. Left: linear scale. Right: 
logarithmic scale. 




Figure 23: Cross section for Z' production in pp collisions at 14 TeV. Left: linear scale. Right: 
logarithmic scale. 
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Table 21: LO Z' production cross section in pb at the LHC for pp collisions at y/s = 1 TeV 
for the various models in Table [Hand varying from 1 to 5 TeV, with steps of 500 GeV. The 
CTEQ6L LO parton distribution functions are employed. 









<*(4) 


a{Z[) 


o (4) 




a (4sM) 


1.0 


0.57 


0.38 


0.41 


0.70 


0.88 


0.89 


1.59 


1.5 


6.3 x 10~ 2 


4.2 x 10~ 2 


4.5 x 10~ 2 


7.0 x 10~ 2 


8.6 x 10~ 2 


8.7 x 10~ 2 


0.16 


2.0 


7.7 x 10~ 3 


6.1 x 10~ 3 


6.4 x 10~ 3 


8.8 x 10~ 3 


1.0 x 10~ 2 


1.0 x 10~ 2 


2.1 x 10~ 2 


2.5 


1.0 x 10~ 3 


9.6 x 10~ 4 


9.8 x 10~ 4 


1.2xl0 -3 


1.3 x 10~ 3 


1.3 x 10~ 3 


2.9 x 10~ 2 


3.0 


1.5 x 10~ 4 


1.4 x 10~ 4 


1.4xl0 -4 


1.5 x 10~ 4 


1.5 x 10~ 4 


1.5 x 10~ 4 


3.9 x 10~ 4 


3.5 


1.7 x 10~ 3 


1.9 x 10~ 3 


1.8 x 10~ 3 


1.7 x 10~ 3 


1.5 x 10~ 3 


1.5 x 10~ 3 


4.7 x 10~ 3 


4.0 


1.7 x 10-" 


1.9 x 10~ b 


1.8 x 10" 6 


1.5 x 10~ b 


1.4 x 10-" 


1.3 x 10~ b 


4.4 x 10~ b 


4.5 


1.1 x 1(T ' 


1.3 x 10 -7 


1.2 x 10 -7 


1.0 x 10"'' 


9.2 x 10~ s 


9.1 x 10~ 8 


3.0 x 10" 7 


5.0 


5.5 x 10~ y 


6.0 x 10~ y 


5.9 x 10~ y 


5.1 x 10~ y 


4.6 x 10~ y 


4.5 x 10~ y 


1.4 x 10~ 8 



Table 22: As in Table I2T1 but at the centre-of-mass energy y/s = 8 TeV. 



m z < 




<y{z' w ) 


<*(4) 


o(Z[) 


o (Z's) 


o(Z> x ) 


a (4sM) 


1.0 


0.83 


0.54 


1.01 


1.28 


1.28 


1.30 


2.30 


1.5 


0.10 


6.9 x 10~ 2 


7.4 x 10~ 2 


0.12 


0.15 


0.15 


0.27 


2.0 


1.6 x 10~ 2 


1.2 x 10~ 2 


1.2 x 10~ 2 


1.8 x 10~ 2 


2.2 x 10~ 2 


2.2 x 10~ 2 


4.3 x 10~ 2 


2.5 


2.8 x 10~ 3 


2.3 x 10~ 3 


2.4 x 10~ 3 


3.1 x 10~ 3 


3.5 x 10~ 3 


3.5 x 10~ 3 


7.4 x 10~ 3 


3.0 


4.9 x 10~ 4 


4.5 x 10~ 4 


4.6 x 10~ 4 


5.2 x 10~ 4 


5.5 x 10~ 4 


5.6 x 10~ 4 


1.3 x 10~ 3 


3.5 


8.4 x 10~ 3 


8.5 x 10~ 3 


8.5 x 10~ 3 


8.4 x 10~ 3 


8.3 x 10~ 3 


8.3 x 10~ 3 


2.2 x 10~ 4 


4.0 


1.3 x 10~ 3 


1.4 x 10~ 3 


1.4 xlO - * 


1.2xl0~ s 


1.1 x 10~ 3 


1.1 x 10~ 3 


3.4 x 10~ 3 


4.5 


1.7 x 10~ b 


1.9 x 10~ b 


1.8 x 10~ b 


1.6 x 10~ b 


1.4 x 10" b 


1.4 x 10~ b 


4.4 x 10~ b 


5.0 


1.7 x 10~ 7 


1.9 x lO" 7 


1.9 xlO -7 


1.6 xlO -7 


1.4 x 10~ 7 


1.4 xlO -7 


4.5 x 10- ' 



Table 23: As in Tables I2T1 and l22l but at the centre-of-mass energy y/s = 14 TeV. 





<*(4) 


a{Z' w ) 


<*(4) 


a{Z[) 


<*(4) 


a(4) 


a (4sM) 


1.0 


2.87 


1.83 


2.00 


3.56 


4.53 


4.59 


7.98 


1.5 


0.52 


0.34 


0.37 


0.64 


0.82 


0.83 


1.43 


2.0 


0.13 


9.0 x 10~ 2 


9.7 x 10~ 2 


0.16 


0.20 


0.21 


0.36 


2.5 


4.0 x 10~ 2 


2.8 x 10~ 2 


3.0 x 10~ 2 


4.9 x 10~ 2 


6.1 x 10~ 2 


6.2 x 10~ 2 


0.11 


3.0 


1.3 x 10~ 2 


9.7 x 10~ 3 


1.0 x 10~ 2 


1.6 x 10~ 2 


2.0 x 10~ 2 


2.0 x 10~ 2 


3.6 x 10~ 2 


3.5 


4.8 x 10~ 3 


3.6 x 10~ 3 


3.8 x 10~ 3 


5.5 x 10~ 3 


6.6 x 10~ 3 


6.7 x 10~ 3 


1.3 x 10~ 2 


4.0 


1.7 x 10~ 3 


1.4 x 10~ 3 


1.4x 10~ 3 


2.0 x 10~ 3 


2.3 x 10~ 3 


2.3 x 10~ 3 


4.6 x 10~ 3 


4.5 


6.4 x 10~ 4 


5.4 x 10~ 4 


5.6 x 10~ 4 


7.1 x 10~ 4 


8.0 x 10~ 4 


8.1 x 10~ 4 


1.7 x 10~ 3 


5.0 


2.4 x 10~ 4 


2.1 x 10~ 4 


2.2 x 10~ 4 


2.5 x 10~ 4 


2.8 x 10~ 4 


2.8 x 10~ 4 


6.2 x 10~ 4 
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foreseen number of events in Z ! decays is simply given by the product of integrated luminosity, 
production cross section and relevant branching ratio. 

The expected event rates in the two considered scenarios are summarized in Tables 1241 and 
l25l for m z /=1.5 and 2 TeV and setting the slepton mass m~ to the values which in Tables [8]- 
l20l maximize the slepton rate. We discarded the Z' x model as it does not yield a sfermion 
spectrum after the addition of the D-term to squark and slepton masses. As discussed in Section 
[5] leptonic final states in supersymmetric events can be yielded by direct decays Z' — > £ + £~ 
(Fig.[T3l) or by a cascade originated from primary decays into sneutrino, chargino or neutralino 
pairs (see Figs. [14] [15] and [161). By adding up such rates, one obtains the so-called cascade 
branching ratio: 

BR casc = BRvv* + BR^ +r + BR |0| o . (34) 

In Tables [24] and [25] N s i ep and N casc are the number of events with a Z' decaying into a 
primary charged- slepton pairs or into a supersymmetric cascade, respectively. In both luminos- 
ity (energy) regimes, due to the large cross section, the Sequential Standard Model is the one 
yielding the highest production of supersymmetric particles in Z' decays, up to ^(10 4 )-<^(10 5 ) 
for cascade events at y/s = 14 TeV and /Jzfdf=100 fb 1 and a Z' mass m z i = 1.5 TeV. As 
discussed in Sections [5TT1 and [541 in the Z^ and Z[ models direct decays into charged sleptons 
are prevented, but sneutrino, neutralino and chargino productions are accessible, with expected 
number of events from 50 to ^(10 4 ) according to m Z ', energy and integrated luminosity. As 
for the Z' N model, in the high-luminosity phase, a few hundreds of direct sleptons and up to 
10 4 cascade particles can be produced for m z , = 1 .5 TeV. For /jSf d?=20 fb 1 and y/s = 8 TeV, 
direct slepton decays are negligible, but about 400 and 70 cascade events can be expected for 
a Z' mass of 1.5 and 2 TeV, respectively. The Z' s boson leads to many cascade particles in the 
high-luminosity regime, between 10 3 and 10 4 , and a few tenths of direct leptons. For the lower- 
luminosity case, there are no directly produced charged sleptons, whereas the cascade sparticles 
are about 30 (m Z '=1.5 TeV) and 46 (m z / = 2 TeV). 

Before concluding this subsection, we point out that, although the numbers in Tables [24] 
and [25] encourage optimistic predictions on Z' decays into sparticles, especially in the high- 
luminosity phase, before drawing a conclusive statement on this issue, it will be necessary car- 
rying out a careful study accounting for detector acceptance and resolution, triggering efficiency 
and cuts on final-state jets and leptons. Hence, the results presented in this paper should be seen 
a first step towards a more thorough investigation, which requires, above all, the implementa- 
tion of the models herein discussed into a Monte Carlo event generator. In this perspective, 
one should compare the Monte Carlo predictions with the experimental data following, e.g., the 
approach proposed in [[441 or investigating the observables suggested in [|45ll to search for new 
physics in Drell-Yan like events mediated by a new heavy resonance. The same analysis should 
be also performed for the Standard Model backgrounds: as discussed in the Introduction, we 
shall defer the implementation of the modelling for Z' production and decays, as well as the 
comparison with the simulation of the backgrounds, to future work. 
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Table 24: Number of supersymmetric particles at the LHC, for Z' production U(l)' models 
and in the Sequential Standard Model at y/s=& TeV and / =Sf d?=20 fb _1 , as a function of m z t 
expressed in TeV. 



Model 




N 

1,1 case 


N s lep 


Z' 


1.5 


523 




Z' 


2 


55 




z' 


1.5 


599 


36 


z' 


2 


73 


4 


7' 

N 


1.5 


400 


17 


7' 

N 


2 


70 


3 


7' 


1.5 


317 




7' 


2 


50 




7' 
^s 


1.5 


30 




7' 
^s 


2 


46 




Z 

^SSM 


1.5 


2968 


95 


7' 

A SSM 


2 


462 


14 



Table 25: As in Table[2l but for ^ = 14 TeV and fSf d?=100 fb _1 . 



Model 


m z / 


N 

^casc 


Nslep 


Z' 


1.5 


13650 




Z' 


2.0 


2344 




z 


1.5 


10241 


622 


z' 


2.0 


2784 


162 


7' 

N 


1.5 


9979 


414 


7' 

N 


2.0 


2705 


104 




1.5 


8507 






2.0 


2230 




4 


1.5 


8242 


65 


7' 
A s 


2.0 


2146 


16 


7' 

A SSM 


1.5 


775715 


24774 


7' 

A SSM 


2 


19570 


606 
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7 Conclusions 



In this paper, we discussed production and decay of new neutral Z' bosons, according to new 
physics models based on a U(l)' gauge group and in the Sequential Standard Model. Unlike 
most analyses undertaken so far, based on SM decays, we also included Z' supersymmetric de- 
cay modes, as predicted by the Minimal Supersymmetric Standard Model: in this perspective, 
the current Z' mass limits may have to be revisited. Extending the MSSM with the U(l)' sym- 
metry implies new features, such as an extra scalar neutral Higgs boson, two novel neutralinos 
and a modification of the sfermion masses due to an additional contribution to the so-called 
D-term. The particle mass spectra were studied in terms of the parameters characterizing the 
U(l)' group and the MSSM; in particular, we discarded scenarios wherein, for fixed values of 
the sfermion soft mass, squarks or sleptons are not physical after the addition of the D-term. 
The same study has been performed for the purpose of the Z' partial widths and branching ra- 
tios, paying special attention to final states with charged leptons and missing energy. In fact, 
these configurations are favourable for an experimental detection at hadron colliders and can 
be yielded by intermediate charged sleptons or a supersymmetric cascade through neutralinos, 
chargino and sneutrinos. The branching ratios of these Z' decays have been investigated in all 
the models, as a function of the slepton mass. 

We finally computed the Z' production LO cross section in all scenarios and gave an estimate 
of the number of supersymmetric events in Z' decays, in the narrow-width approximation and 
for few values of centre-of-mass energy and integrated luminosity. The outcome of this study 
is that, for some models and parametrizations, one can even have up to 10 4 -10 5 events with 
sparticle production in Z' decays. As an additional remark, we wish to point out that the Z' — > 
£ + £~ decay presents two interesting features. First, the Z' mass will set an additional constrain 
on the slepton invariant mass; second, it allows the exploration of comers of the phase space 
which would be instead unaccessible through other processes, e.g. Drell-Yan like events. 

In summary, we consider our investigation a useful starting point to study Z' production and 
decay beyond the Standard Model, such as within supersymmetric theories, drawing guidelines 
for future experimental analyses. In future perspective, it will be very interesting performing 
a study including parton showers, finite-width and hadronization corrections, as well as ex- 
perimental effects, like the detector simulation and the acceptance cuts. In this way, one will 
eventually be able to draw a statement on the Z' mass limits within supersymmetry. To reach 
these objectives, the models for Z' production and decay, examined throughout this paper, will 
have to be implemented in Monte Carlo programs, such as HERWIG or PYTHIA, and the super- 
symmetry signals compared with the Standard Model backgrounds simulated, e.g., by means of 
the ALPGEN code 11461 . In the framework of an event generator, it will also be possible, in the 
same manner as the experimental analyses do, rescaling the total cross section in such a way to 
include higher-order QCD corrections. For this purpose, the use of the FEWZ code [|47l . which 
simulates vector boson production at hadron colliders at NNLO, with fully exclusive final states, 
is planned. Other possible extensions of our analysis consist in investigating more thoroughly 
the unconventional assignment of the SM and exotic fields to the SU(10) representations, as 
well as scenarios wherein the exotic leptons (sleptons) and quarks (squarks), predicted by the 
grand-unified group Eg, but discarded in the present work, are lighter than the Z' and therefore 
capable of contributing to its decay width. This is in progress. 
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A 2J decay rates into standard and supersymmetric particles 



The Lagrangian term describing the interaction of the Z' with fermions is given by: 

J? f = g'ff(v f -a f y 5 )fZ^ 

with 

Setting Q'(f R ) = —Q'(f[), the vector and axial-vector couplings read: 

v f = ~ [Q'(f L ) + Q'(f R )] ,a f =\ [Q'(f L ) - Q'(f R )] , 



where the U(l)' charges of left- and right-handed fermions can be obtained by using Eq. £[2 
and Table [2] In terms of the mixing angle 9, such couplings read: 



(35) 
(36) 

(37) 



r t - - (Q' v (fL)+Q' v (fR))cose-(Q' x (f L ) + Q' x (f R ))sme 
a f = ^ [(Q' v (f l) - Q' v ( f r)) cos 9 -(Q' z (f L )-^ x (f R )) sin 9 
One can thus write the Z' width into fermion pairs as: 



(38) 



r(z'^ff) = c f 



J2 



\2% 



m z , 



fj2~ \ f TY\^ 

1 +2^- I +a} 1-4-f 



m- 



1-4 



/// 



1/2 



f_ 



(39) 



"z 1 J \ " l z' 

where the colour factor is Cf = 3 for quarks and Cf = 1 for leptons. With the charges listed in 
Tableland employing Eq. (|38l , one can show that, in the Z[ model, namely 9 = arccos a/5/8 — 
7t/2, the vector and vector-axial couplings of the Z' with up-type quarks vanish, i.e. v u = a u = 0. 
In fact, when discussing Z[ phenomenology at the Representative Point (Section 4), it was 
pointed out that its branching ratio into uu pairs is null. 

Likewise, the interaction Lagrangian of the sfermions with the Z' reads: 



X f = g'{v f ±a f ) [fUd^x) - (d^ R )f L>R }Z'». 
The width into left- or right-handed sfermions is given by: 

Jl ( mY V2 

r(Zf -+ f LM fi R ) = Cff^m z , (v f ± a f f f 1 - 4-f 



(40) 



(41) 
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where the ± sign refers to left- and right-handed sfermions, respectively. Eq. (I4TT) is expressed in 
terms of weak eigenstates fiR, its generalization to the mass eigenstates /i 2 is straightforward 
and discussed in lfT2ll . However, for the parametrizations used throughout this paper, sfermion 
mixing is always negligible and Eq. (HTI) can be safely used even to calculate the branching 
ratios into f\fl and /2/I final states. 

From Eq. (RTb one can learn that the Z' rate into left- and right-handed sfermions vanishes 
for Vf = —cif and vy = at, respectively. In fact, for vy = a/, according to Eq. (|35l) . the Z' only 
couples to left-handed fermions and therefore in the MSSM, in absence of left-right mixing, 
there is no coupling with right-handed sfermions. Likewise, for Vf = —cif, the Z' only couples 
to right-handed fermions and sfermions and the rate into fif[ pairs is null. For example, in the 

model, it is v v = a v , whereas, in the Z[ model, V| = a~ t Therefore, as remarked in Sections 
5.3 and 5.4, the — > V2V2 an d Z[ ~~ ^ ^2^2 are suppressed, although they are kinematically 
permitted at the Reference Point. 

As for the Higgs sector, defining Q\, Q' 2 and Q' 3 the U(l)' charges as in Eq. (fT5l and /3 = 
arctan(v2/vi), one can obtain the Z' rate for decays into charged-Higgs pairs [48] 



12 ( 2 \ 3 / 2 

T(Z'^H+H ) = A_ ( 2 ;sin 2 /3 -Q 2 cos 2 /3) 2 m z A\ (42) 



and associated production of a W boson with a charged Higgs |^ 



T(Z'^W ± H^) = ^( e ; + e ' 2 ) 2 mz ,sin 2 /3cos 2 /3 

487T 



2 ^2 



| 5nfy-mjj± (m 2 w -m 2 H± 

2 + 4 



1 2 <+ ffl2 j± , ( m w- m2 H ±)\ (43) 
m 2 , m z 

As the Z' has no direct coupling with W's, the decay into pairs occurs by means of the 

Z-Z' mixture. For small values of the Z-Z' mixing angle, this width reads Il48l : 

T(2f = ^ (e'jcos 2 ^ -2 2 sin 2 /3) 2 m z ,. (44) 

In order to obtain the widths into Z-Higgs pairs, i.e. Zh, ZH or ZH' final states, or into scalar- 
pseudoscalar neutral-Higgs pairs, such as hA, HA or H'A, one first needs to diagonalize the 
neutral Higgs mass matrix (see llT2lO . The Z-Higgs rate can be written in compact form as: 



e '2 

T{Z'^Zhi) = ^(Q[cosPU u -Q 2 sinl3U2i) 2 m z , 



2 \2' 



, ^ m z- m hi ( m w- m h<) 



m z , m z , 



x ^.^tH-^ 1 (45) 

m z , m z , 



9 ^Eq. d43l > corrects a typing mistake present in Ref. [12], wherein the decay width Z' — > W ± H^ : is instead 4 
times smaller than in ( l43l >. 
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where i = 1,2,3 for final states Zh, ZH and ZH', respectively, and Uij is the matrix which 
diagonalizes the Higgs mass matrix in the (h H H') basis. Likewise, using the same notation 
as in Eq. (|45T ), the scalar-pseudoscalar Higgs width reads: 



T(Z' -> h t A) 



x 



12 2 



487T/V 2 



{v 3 Q[ smpUu + v 3 & 2 cospU2i + vQ' 3 sin/3 cos/3t/ 3; ) 2 m z , 



1-2 



m l + m A 



+ 



m 



2\2' 



nv 



3/2 



(46) 



In Eq. (l46l) . following [fT2|. we defined N = y v\v\ + v 2 v 2 + v 2 v 2 and v = J v\ + v\. 

Finally, one can derive the decay widths into gauginos. As for neutralinos, after diagonaliz- 
ing the mass matrix (|2TT ). the interaction Lagrangian reads: 



(47) 



hj 



where gjj is a generalized coupling depending on the diagonalizing-matrix elements and has 
been calculated numerically. The partial rate into neutralino pairs (%f%j) with masses m, and 
nij is thus given by: 



r(z' fix)) 



g 



\2% 



mj+m: 



m 



2\2 



2m 7 , 



2m%, 



mjtnj 

2 

mi, 



x 



\ 



1- 



(mi + mj] 



iik. 



mj) 



(48) 



Finally, the Lagrangian term corresponding to the coupling of the Z' with charginos is given by: 



(49) 



hi 



The generalized vector and vector-axial couplings can be expressed in terms of 0±, the angles 
of the unitary transformation diagonalizing the chargino mass matrix [18], and the Higgs U(l)' 
charges as follows [fT2|: 



vii = ! 


2 / 1 sin 2 ^-e 2 sin 2 + , 




an = ! 


2isin 2 0_ + 2 2 sin 2 + , 




V\2 = V21 = ! 


2j sin 2 0_ cos 0_ — 5<2 2 s i n 0-+ 


_ -fcos (j) 


«12 = «21 = 


2 1 sin 2 0_ cos + + 5<2 2 sm 0-+ 


- +COS 


V22 = 


2iCOS 2 0_-2 2 cos 2 + , 




«22 = 


2iCos 2 0_ + 2 2 cos 2 0+. 





(50) 



In the above equations, 8 = sgn(my±)sgn(m^±) . The analytical expressions for <j>± can be found 



*2 
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in [fT8l and are not reported here for brevity. The rate into chargino pairs is finally given by: 



J2 



487T 



m z , <J {vfj + ajj) 



mf + mj (mf — mj 
1 J 1 j 



2\2' 



2mL 



x 



\ 



_ (mj + nij 



■ \2 



m z> 



m z> 



2 minij 
mf, 



(51) 
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